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SYNOPSIS 

The continuously growing demand for electric power requires the transmission of large 
amounts of po'v\er over long distances An economically attractive solution to increase 
power transfer through long transmission lines, without building new parallel circuits, is 
to install the series capacitors It is known that series capacitor compensation benefits 
power systems in many wa’vs, such as enhancing transient stability limits, increasing power 
transfer capability, etc It is also known that fixed senes compensation may cause sub 
synchronous resonance (SSR) m power systems, which can lead to the damage to machine 
shaft Since the discovery m 1970 that SSR was the mam cause of the shaft failures at the 
Mohave generating station (USA) extensive research and development efforts have been 
devoted to the development of effective SSR mitigation measures 

The definition of SSR is given by IEEE [1] Subsvnchronous resonance is an lelectnc 
power syster^condition where the electric network exchanges energy with a turbine genera 
tor at one or more of the natural frequencies of the combined system below the synchronous 
frequency of the system 

When a series compensated power system is perturbed, its equilibrium state is dis 
turbed, giving rise to interchange of energies between mechanical and electrical systems 
mutually coupled through the rotor of the synchronous generator As energies are inter 
changed, their frequencies of oscillations are the natural frequencies of their respective 
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systems The natural frequencies of oscillation of each system form their respective modes 

of oscillation An oscillation of generator rotor, due to disturbances, produce super syn 

chronous and sub synchronous components of currents (or voltages) These currents induce 

subsynchronous and supersynchronous torques There are many ways m which the system 

and the generator may interact with subsynchronous effects Couple of these interactions 

are (1) Torsional interaction and (2) Transient torque effect The turbine generator shaft 

system responds to disturbances with oscillations at its natural frequencies Oscillations of 

the geneiator rotor at its natural frequency result m modulation of generator voltage The 

subsvnchronous frequency voltage is the complement of the natural frequency of the turbine 

generator shaft system When this subsynchronous frequency is close to system natural 

frequency, the resulting armature currents produces a subsynchronous torque which re 

f 

mforces the aforementioned generator rotor oscillations This can result in sustained or 
growing oscillations This phenomenon is referred to as torsional interaction 

System disturbances impose electromagnetic torques on generator rotor, subjecting 
shaft segments to torsional stresses Following a significant system disturbance in series 
capacitor compensated power system, the resulting electromagnetic torque oscillations at a 
frequency complement to the electrical system natural frequency If this frequence coincides 
with any of the natural modes of the turbine generator shaft, there can be large peak 
torques This effect is referred as transient torque effect or shaft torque amplification 
Power electronic systems are finding increasing applications in power systems for both 
transient state and steady state operation requirements [2] STATCOM is a second gener 
ation FACTS device The STATCOM consists of GTO based voltage source inverter (VSI) 
with a dc capacitor coupling transformer and control circuit The VSI is connected to the 
AC bus through coupling transformer The variation in STATCOM output voltage can be 
achieved through the control of the firing angle [3] STATCOM with voltage controller for 
SSR mitigation is proposed in the literature The voltage controller along with the reactive 
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current controller may be considered for improving the performance of the ST\TCOM 

Controlled series compensation (CSC) is a type of Flexible A.C transmission systems 
(FACTS) controller which utilizes thyristor valves to control the degree of series compensa 
tion [4] CSC consists of tvo components, one element is a mechanically switched portion 
the second portion being a Thyristor Controlled series capacitor (TCSC) The TCSC con 
sists of a fixed capacitor in parallel with a thyristor controlled reactor(TCR) The TCR 
reactance is varied by varying the firing angle of the thyristors, which m turn varies the 
effective capacitive reactance Discrete time and continuous time TCSC models were pro 
posed in the literature to study the SSR phenomenon and to estimate the damping of 
torsional modes in open loop The TCSC control strategy can be utilized to study the 
impact of TCSC on torsional modes 

It IS mentioned previously that senes compensated power system gives raise to subsyn- 
chronous phenomenon Vanous methods of counteracting subsynchronous resonance phe 
nomenon have been considered which are great diversity in concept and methods Most of 
these concepts were tested with IEEE first bench mark model One of these methods is to 
modulate either real power or reactive power injected into the generator bus of the affected 
steam turbine The modulation of real power and reactive power is one of the effective 
means of mitigating subs nchronous resonance and this can be achieved by FA.CTS con- 
trollers The real power modulation can be achieved by the modulation of TCSC reactance 
wheieas reactive power modulation can be obtained by STATCOM 

The tmsional interaction may be considered as small disturbance phenomenon whereas 
transient torque problem is due to the large disturbance The torsional interaction analysis 
IS usually carried out utilizmg linearized models The study of transient torque problem 
requires the detailed modeling of all system nonlinearities Digital simulation is well suited 
for analyzing transient torque effect [1] 

Subsynchronous resonance control through FACTS controllers is gainmg importance 
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Based on the above considerations, this thesis focussed the attention on control of subsyn 
chronous resonance which includes torsional interaction and transient torque effects with 
FACTS controllers such as STATCOM and TCSC 
The objectives and scope of thesis are 

1 To study the subsynchronous transient oscillation damping with the ST ATCOM 

2 To propose a new digital control scheme for TCSC control to damp subsvnchronous 
resonance oscillations 

3 To study the control mteractions between TCSC and power system stabilizer 

An outline of the work reported in the thesis is given below 

The first chapter gives an introduction to the various aspects of the problem presented 
m the thesis and reviews bnefly the previously published literature 

In Chapter 2 linearized model of generator system with detailed representation of 
stater, rotor and mechanical system along with AC network is presented The study 
system in this Chapter and subsequent Chapters is the celebrated IEEE First bench mark 
model [1] Eigenvalue anal} sis is carried out with the linearized overall system model 
The results obtained through eigemalue analysis is validated with time domain simulation 
and Fourier analysis The effectiveness of NGH damping scheme [5,6] in mitigating the 
subsynchronous torque oscillations is studied through the detailed time domain simulation 
In Chapter 3, linearized model of 12 pulse STATCOM along with voltage and reactive 
controller is presented [7] The control parameters are obtained through eigenvalue analysis 
of the linearized model The results obtained through eigenvalue analysis are validated with 
time domain simulation for the study of transient torque oscillations Comparative study 
IS carried out with NGH damping scheme in mitigating subsynchronous torque oscillations 
SSR analysis with TCSC m continuous time linearized domain is presented m Chapter 
4 Constant angle control of TCSC is investigated The constant angle control is based on 
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the philosophy of maintaining the voltage drop across the compensated transmission line 
constant Eigenvalue anal} sis is earned out to predict the stability of the system and the 
results are validated with time domain simulation 

In Chapter 5 a no'vel digital control scheme is proposed for TCSC control For this 
purpose a discrete time domam TCSC model is developed The controller is designed using 
state feedback approach The discrete time model along with the controller is suitably 
interfaced with the generator system model to obtain a homogeneous state space equation 
The stability of the system is analyzed through eigenvalue analysis and the results are 
validated with time domain simulation 

The effect of excitation system that includes power system stabilizer (PSS) on tor 
sional modes, in the presence of TCSC controller is examined in Chapter 6 The control 
interaction study is carried out through eigenvalue analysis 

This thesis concludes in Chapter 7 outlining the conclusions drawn from the thesis and 
suggest some future scope of woik 
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Chapter 1 

INTRODUCTION 


One of the important challenges facing the utility industry is the efficient utilization of 
the existing transmission network, in view of the difficulties invohed in adding new trans 
mission capacit\ The increase in power flow over a given transmission network can be 
achieved by compensating the AC network either by senes compensation of part of indue 
tive reactance of the transmission line by series capacitors or bv shunt compensation to 
maintain voltage at appropriate points yi the network bv reactne power compensators 
Series capacitors provide a direct approach to increase transmission capability which 
IS very often the most economical solution Unfortunately however, it can give rise to 
subsynchronous resonance (SSR) by interacting with the turbo-generator Two incidents 
of generator shaft damage at Mohave generating station in the Umted States [1] prompted 
investigation into the nature of the problem of subsynchronous resonance and its solution 
The slow nature of series capacitors control through mechanical swit chin g restricted their 
usefulness which implies that faster dynamic controls are required for better utilization 
of AC transmission network The advent of high voltage and high power thyristor valves 
and digital controllers in HVDC transmission has demonstrated the viability of fast con- 
trollers for power transmission subsequently this has led to the concept of Flexible AG 
transmission systems (FACTS), which uses advanced solid state controllers to overcome 
the problem of underutilization of AC transmission networks 
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FACTS includes host of controllers namely static var compensator (SVC), thyristor 
controlled series capacitor (TCSC), static compensator (STATCOM), thyristor controlled 
braking resistor (TCBR), thyristor controlled phase angle regulator (TCPR) unified power 
flow controller (LPFC), and solid state series compensator (SSSC) In this thesis an attempt 
has been made to apply some of the FACTS controllers to control the subsynchronous 
resonance caused bv series capacitors Subsynchronous resonance phenomenon and the 
various FACTS controllers are discussed in the following sections 

1 1 Subsynchronous Resonance 

Subsynchronous resonance is a major concern for the stability of turbine generators con 
nected to transmission systems that employ series capacitors A disturbance in the power 
system can cause stimulation of turbine-generator natural torsional modes When the 
generator is connected to a senes capacitor compensated system these oscillations can be 
amplified and sustained due to interaction between the electrical power system and the 
torsional mechanical system SSR has been defined by the IEEE SSR Task Force [2] as 

Subsynchronous resonance is an electrical power system condition where the electrical net 
work exchanges energy with the turbine generator at one or more of the natural frequencies 
of the combined system below the synchronous frequency of the system 

During the incidents of generator shaft damage at Mohave [1], it was found that the 
frequency of one of the torsional modes was close to the complementary frequency of 
subsynchronous currents present in the electrical system This resulted in large torque in 
the shaft section bet veen the generator and exciter which subsequently damaged the shaft 

There are two pr iblems associated with SSR [3] One is steady state SSR which may be 
due to (a) induction generator effect and (b) torsional interaction The second problem is a 
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transient SSR phenomena and termed as transient torque or torque amplification problem 
To understand these aspects, consider a synchronous generator connected to an infinite 
bus through a series compensated line, as shown in Fig 1 1 The series resonant frequency 
of the electrical network is given by 

X +Xl 

where 

X" subtransient reactance of generator 
X[^ reactance of line and transformer 
Xc reactance of series capacitor 
fo nominal frequenc\ 

It IS to be noted that since Xc < Xl, fe < fo 



Rl X l 

n/VV 1 ^ 

Generator 


Figure 1 1 Series compensated power system 



Infinite Bus 


The electrical disturbances that occur in the system cause transient currents which die 
away with a damping rate depending on the value of series resistance Ri at frequency 
/e These subsynchronous currents of frequency fe fiiow in the armature and result m 
subsynchronous torques and currents in the rotor circuit at frequency fr given by 

= ( 12 ) 
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The subsynchronous rotor currents and mechanical oscillations produce subsynchronous 
voltage components which sustain the currents The frequency fr is defined as the comple 
mentary frequency of fe 

111 Induction Generator Effect 

The induction generator effect is due to the electrical system dynamics The rotor resistance 
to subsynchronous currents viewed from the armature terminals is negam e A.t the system 
resonant frequencv fe, if this negative resistance exceeds the positive resistance of network, 
self excitation takes place However this problem can be tackled bv suitable design of 
amortisseur circmts [4] 

112 Torsional Interaction 

The Generator rotor oscillations at a torsional mode frequency, fm induce armature voltage 
components at frequencies (fem) given by 

/em = /o ± fm (1 3) 

When fem IS close to fe (electrical resoqajit frequency), the subs-v nchronous torques pro- 
duced by the subs\ nchronous voltage component can be sustained [4] This mterplay 
between electrical and mechanical systems is termed as torsional mteraction 

113 Torque amplification Problem 

System disturbances impose electromagnetic torques on generator rotors, subjecting shaft 
segments to torsional stresses Following a significant disturbance in a series capacitor 
compensated system, the resulting electromagnetic torque oscillates at frequency of fo~fe 
If this frequency is near to any natural frequency of shaft section, the resultmg shaft torques 
could be much larger than those produced by a three phase fault in a system without series 
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capacitors This is caused by resonance between the electrical system and mechanical 
system natural frequencies These effects are referred to as shaft torque amplification [5] 

114 Countermeasures for SSR 

The control of SSR during transients in senes compensated systems using static filters 
bypassing of senes capacitors, excitation control etc , have been suggested in literature 
[6-8] A summary of the various solutions for mitigating the effect of SSR is presented in [9] 
Farmer et al [8] Fouad et al [7] and Li Wang [10] suggested the injection of a properly 
phased sinusoidal signal from the rotor motion into the excitation system Edris [11] 
proposed a method of counteracting SSR based on the decoupling of machine network 
interactions The decoupling is achieved by introducing senes capacitor compensation 
scheme which behaves as a set of three phase asymmetrical capacitors at subsynchronous 
frequencies Ira\mii et al [12] proposed to modulate generator active power by injecting 
a quadrature phase voltage with thyristor controlled phase shifter Hammad et al [13] 
suggested thyristor phase controlled reactor utilizing its mam regulator and an auxiliary 
speed signal to damp all torsional modes The modulation of real power and reactive power 
is one of the effectne means of mitigating the subsynchronous resonance 

Hingoram et al [14, 15] have proposed the NGH damping scheme to reduce transient 
machine oscillations during system disturbances The scheme involves a hnear resistor in 
series with a pair of anti parallel thyristors connected across the capacitor When a zero 
voltage crossing point of the capacitor voltage is detected, the succeeding half cycle period 
is timed If and when the half cycle exceeds the set time, the corresponding thyristor is 
fired to discharge the capacitor through the resistor and bring about its current zero sooner 
than otherwise The thyristor stops conducting when the current through it reaches zero 
The modulation of real power and reactive power is one of the e fective means of 
mitigating the subsynchronous resonance and this can be achie\ed by FACTS conti oilers 
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The real power modulation can be achieved by the modulation of TCSC reactance whereas 
the reactive power modulation can be obtained by STATCOM A review of various FACTS 
controllers is given below 

1 2 Review of FACTS Controllers 

The power flow in an AG line is dependent on three important parameters namely voltage 
magnitude impedance and phase angle FACTS devices control one or more of the param- 
eters for improving the transmission capability Owing to the accuracy and fast response, 
FACTS devices are used to improve the dynamic performance, under steady state as well 
as transient state 

12 1 Static Var Compensator (SVC) 

SVC IS the first generation FACTS device that is alreadv in operation at various places 
around the world It uses a fixed capacitor (FC) or thvristor switched capacitor (TSC) 
m conjunction with a thyristor controlled reactor (TCR) SVC controls the magnitudes of 
the voltages at selected buses in the transmission network For example locating a SVC at 
the midpoint of a transmission line can increase power flow in the line The mam feature 
of SVC IS the voltage control by means of the reactive power compensation (obtained by 
varying the firing angle of the thyristors) In addition, SVC may improve transient stability 
by dynamically supporting the voltage at key points and steady state stability by helping 
to increase the swung mode damping SVG with proper control signals can be used to 
damp out electro mechanical oscillations Damping of power system oscillations play an 
important role not only in mcreasmg the transmission capability of the line, but also for 
stabilization of S} stem after critical faults in weakly coupled systems 
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12 2 Controlled Series Compensation (CSC) 

CSC provides a wav to increase the power flow on selected lines without causing the 
problems previoush associated with series capacitors It is a second generation FACTS 
controller which controls the impedance between the two buses to which it is connected 
and hence the effectue line reactance CSC consists of conventional series capacitors which 
are mechanicallv switched and a thyristor controlled series capacitor (TCSC) A single 
TCSC module consists of a series capacitor and a parallel path with a pair of anti parallel 
connected thyristor \alves and a series inductor Also in parallel is a metal oxide varistor 
(MOV) for over toltage protection A complete compensation system mav be made up of 
several of these modules m series 

The TCSC module has three basic operating modes (i) Bypassed mode (ii) Inserted 
with thyristor valve blocked (iii) Inserted with vernier control In the bvpassed mode 
thyristors are gated for full conduction (180°) and the current flow in the reactor is contin 
uous and sinusoidal The net reactance is slightly inductive because the susceptance of the 
reactor is larger than that of the capacitor This mode is used mainly for protecting the 
capacitor against the over voltages (during transient ov er currents in the line) This mode 
IS also termed as rhvTistor reactor mode In the inserted mode no current flows through 
the valves with the blocking of gate pulses In this case TCSC reactance is same as that 
of the fixed capacitor and there is no difference in the performance of TCSC in this mode 
with that of a fitxed capacitor This mode is also termed as waiting mode In the inserted 
mode with vernier control, the thyristor valves are gated such that they conduct for a part 
of a cycle (amm < a < 180°) The effective TCSC reactance (in the capacitive region) 
increases as a is reduced below 180° and the TCSC reactance is maximum when a = Orfun 
which may be three times of fixed capacitive reactance This angle (Q!,n,„) is above the 
value of a corresponding to the parallel resonance of thyristor controlled reactor and the 
tixed capacitor (at fundamental frequency) 
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Control of TCSC 

The TCSC control system arrangement is shown in Fig 1 2 [16] It can be observed from 
this figure that the TCSC control functions are partitioned into two le\els i e, common 
level and module level The control commands flow from common level to module level 
while the status information is sent back from each module level The module controller 
executes the ordered change in reactance level The common control lev el receives signals 
of line current and TCSC voltage to generate feedback signals for closed loop control func 
tions It also receives commands from energy management center for settmg power order 
The major control functions are usually power scheduling control power swing damping 
control and transient stability control 

12 3 Static Compensator (STATCOM) 

A STATCOM IS a second generation FACTS controller used for reactive power control and 
incorporates Voltage Source Converters- {VSC) A sue pulse VSC consists of six Gate Turn- 
Off (GTO) thyristors with inverse connected diodes connected as a Graetz bridge with a 
capacitor connected on the DC side The STATCOM is connected to high voltage AC bus 
through a coupling transformer Essentially a STATCOM produces a set of three phase 
output voltages each of which is in phase with and coupled to the respective A.C system 
bus voltage through a small reactance that is generally provided bv the leakage reactance of 
the coupling transformer Ideally as the STATCOM output voltage is m phase with the AC 
system voltage, the reactive power exchange between the AC svstem and the STATCOM 
can be regulated by controlhng the magnitude of the STATCOM output voltage That is, 
if the amplitude of the output voltage is increased above that of the AC system voltage, 
then the current flows through the reactance to the AC system which is tantamount to 
the STATCOM supplying reactive (capacitive) power to the AC system On the other 
hand, if the magnitude of the output voltage of the STATCOM is reduced below that of 
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Figure 1 2 TCSC Control System 
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the AC system voltage then the current flows from the AC SAstem to the STATCOM 
and as a result the STA.TCOM absorbs reactive (inductive) power from the AC bus If the 
magnitude of the output voltage is equal to the AC svstem voltage magnitude, the reactive 
power exchange is zero The converter terminal voltage can be made to lead or lag the AC 
system voltage by controlling the firmg angle [17] Multi pulse converters are usually used 
to reduce harmonics 

12 4 Static Synchronous Series Compensator (SSSC) 

A SSSC IS a solid state voltage source inverter, which injects an almost sinusoidal voltage 
of variable magnitude in series with a transmission hne This has the same structure as 
that of a STA.TCOM except that the coupling transformer is connected m series with the 
transmission line The injected voltage is almost in quadrature with the hne current A 
small part of the injected voltage, which is in phase with the hne current, provides the 
losses m the inverter Most of the injected voltage which is in quadrature with the line 
current emulates inductive or a capacitive reactance m series wurh the transmission hne 
This emulated variable reactance, inserted by the injected voltage source influences the 
electric power flow in the transmission line [18] 

Although subsvnchronous resonance phenomena m senes compensated system has been 
analyzed in detail and countermeasures suggested for overcoming the phenomenon it is 
necessary to analyze the phenomenon in the presence of FACTS devices and their controls, 
as the developments of FACTS devices are quite recent 

1 3 Review of Literature 

13 1 TCSC Modeling 

Most of the TCSC models were det eloped utilizing sampled data approach Ghosh et 
al [19] proposed a disciete time domain TCSC model based on the linearized behavior of 
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the state transition equations that can predict the shift in zero-crossings of the line current 
(or capacitor voltage ) very accurately This model captures the system transient behavior 
accuratelv Rajaraman et al [20] presented a method for computing the eigenvalues of 
the system in discrete time domain The model development m [20] does not present a 
stand alone TCSC model The total system equations were developed with and without 
conductmg of thvTistors through Poincare map which is obtained by integrating the total 
system equations and taking into account the changes in equations when the switching 
occuis This approach requires sampling of full state at even sampling instant which is 
impractical for large power systems 

Othman et al described the TCSC dynamics alone at the switching instants and then 
converted them mto a continuous time model This model assumes that line current is 
forcing function to the TCSC equations and therefore is an mdependent quantity This 
IS the major drawback of the model as change in TCSC firing angle will cause change m 
line current Verghese et al [21] presented a quasi steady state TCSC model based on 
phasor dvnamics approach This model assumes line current is sinusoidal This model 
will be discussed m Chapter 4 Christ et al [22] presented an approach to calculate the 
fundamental frequency impedance of a TCSC This approach will be utihzed later in the 
thesis 

13 2 SSR mitigation with TCSC 

A method for computing the eigemalues and damping of SSR modes based on discrete 
time domain model for the IEEE first benchmark model is presented in [20] The dampings 
associated with subs} nchronous modes of oscillation of a TCSC compensated network are 
computed by linearizing the Poincare map associated with the steady state operation 
The TCSC is cperated viith constant firing angle The analysis is confined to torsional 
modes and impact of TCSC djnamics on network modes is nor mentioned Sampled data 
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model [23] of TCSC along with the higher le\el controllers utilized to stabilize the unstable 
torsional modes of IEEE first bench mark model 

The design installation and basic control modes of the Slatt TCSC is described in 
[24] The Slatt TCSC is comprised of six identical TCSC modules that are connected m 
series A number of field tests were performed to demonstrate the effect of the TCSC m 
mitigating SSR [25] All tests were for radial case and focus on the particular torsional 
mode (mode 4) hich may be destabilized The field test plots show that hile conventional 
series capacitors can cause subsynchronous oscillations the response vhen inserting the 
same capacitor ■v\nth thyristor switch performing vernier control shows no subsvnchronous 
oscillations To e\aluate the ability of TCSC in mitigating SSR caused by a conventional 
series capacitor a small TCSC operating in vernier control was added to a large senes 
capacitor bank ithout the TCSC, mode 4 \\as negatively damped However, the addition 
of TCSC resulted in significant improvement m damping of mode 4 The results of field tests 
were identical to the simulation results obtained using EMTP and TN A. measurements 

The analogue simulator based TCSC studies are presented in [16] The studv examines 
numerous aspects of TCSC performance including response to faults, power swing damping 
and SSR The effectiveness of TCSC was e\ aluated bv measuring the electrical damping 
torque as a function of frequency Comparison of damping torque plots of the system with 
no series capacitor to the svstem with TCSC operating in vernier mode showed that the 
plots are identical i e , the TCSC vernier control resulted in the electrical system look 
like an uncompensated system Comparison of damping torque plots wnth conventional 
series capacitors to that with TCSC operating in vernier mode showed that the TCSC 
significantly reduced the negative electrical damping in the rotor oscillation frequency 
range 

An EMTP study of SSR mitigation of the Slatt TCSC is presented in [26] A dis 
turbance is applied by switching in the series compensation (either fixed compensation 
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or TCSC operating m vernier mode) and the speed deviations of the various masses are 
observed For the case with TCSC torsional modes (1,2 and 4) are significantly smaller 
compared to the case of fixed compensation Frequency domain analysis to compute the 
equivalent TCSC impedance is also presented The TCSC impedance study showed that 
the vernier operation of the TCSC provides equivalent resistive damping 

The Steady state and transient characteristics of the A.SC circuit installed at the 
Kayenta substation are discussed b\ Krause et al [27] The A.SC control mode provides 
controlled firing of the thvnstors to maintain either a fixed capacitive \SC impedance 
(constant reactance control) or to regulate the transmission Ime current The controller 
utilizes a TCSC fundamental frequency impedance verses firing angle look up table for 
converting controller output to thyristor firing angle For reactance control the A.SC con 
trollei determines the proper firing angle to provide the desired reactance In the current 
control mode, a control svstem is used to determine the impedance reqmred to maintain 
the desired line current The electrical damping torque plots (obtained by using a simu 
lator) for cases of a conventional capacitor and an equivalent \SC with reactance control 
showed that while the damping is negative with conventional capacitors the A.SC provides 
positive damping at SSR frequencies 

Bowler et al [28] have described the constant reactance control of TCSC as well as 
modulated reactance control In modulated reactance control, the reactance order is mod 
ulated using an input signal that is s^ nthesized from the voltage and current measurements 
of the capacitor The effectiveness of TCSC in damping power s’svungs and providing higher 
power transfer by permitting higher levels of compensation by mitigating SSR effects is pre- 
sented Damping torque analysis of the system showed that for an unstable SSR condition, 
by adding a small TCSC with constant reactance control reduces the negatne damping 
at torsional mode frequencies Modulated reactance control of TCSC resulted m positive 
damping torque throughout the range of torsional mode frequencies thereby t limmating 
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the destabilizing influence of SSR A root locus plot of the torsional modes as a function of 
controller gam illustrated that the damping of torsional modes is increased with increase m 
controller gam In this study the TCSC has been represented h\ linear transfer functions 
obtained from the time domain simulations 

Heiden et al [29] presented SSR characteristics of different U^pes of series compensation 
schemes The performance of NGH scheme and ad\anced senes compensation scheme 
(ASC) IS analyzed from transient torque and torsional interaction viey point It is shown 
that ASC scheme provides better performance during major disturbances The conclusions 
are drawn based on time domain simulations In the torsional interaction studies it is shown 
that the behavior of ASC is inductive m the subsynchronous frequency range 

13 3 STATCOM and its Characteristics 

The majority of the STATCOM papers published are related to ^oltage control and power 
flow control Edvards et al [30] describe a ±1 MVAR protot\pe Advanced Static VAR 
Generator (ASA G) The purpose of the prototype de\elopment vas to verify performance 
under practical operating conditions This ASVG is based on a 12 pulse voltage source 
inverter The prototype could regulate the line voltage and the reactive power m both the 
leading and lagging directions The mathematical model for the Advanced Static VAR 
Compensator or STATCOM is derived m [31] and frequency domain analysis is used to 
obtain the rele^ant transfer functions for control system synthesis The results are illus- 
trated with measured wa\eforms obtained from a scaled analogue model Schauder et 
al [32] presented the development and installation of a ±100 MVAR STATCOM at the 
Sullivan substation of the TVA power system The operation of the same STATCOM 
IS presented m [17] The effectiveness of voltage control and reactive current control of 
STATCOM IS demonstrated Two Advanced Sta ic Var Compensator (ASVG) configura- 
tions are presented m [33] for power transmission applications The study concludes that 
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the best performance can be achie\ed by ASVC employing PWM switching techniques 
Low frequency oscillation damping with STATCOM is presented in [34] The STA.TCOM 
controller m this case is designed using state feedback approach Time simulations are car 
ried out with 12 pulse STATCOM to verify the interactions benveen the STA.TCOM and 
A.C transmission system The application of STATCOM m damping torsional oscillations 
has not been reported in many papers In [35] an attempt has been made in which the 
the STA.TCOM is connected at the generator terminals It has been sho-nn that in this 
configuration the STATCOM cannot stabilize the torsional modes when operated in the 
voltage control mode The unstable torsional modes are damped with generator speed as 
auxilary signal along with the voltage controller 

1 4 Objectives and scope of the Thesis 

It IS mentioned previously that senes compensated power sjstem gives raise to subsyn 
chronous phenomenon Various methods of counteracting subs^mchronous resonance phe 
nomenon have been considered that have great diversity m concept and methods Most 
of these concepts were tested with IEEE first bench mark model One of these methods 
IS to modulate either real power or reactive power injected into the generator bus of the 
affected steam turbine The use of reactive power modulation is used in [13] to damp sub 
synchronous resonance oscillations The real power and reactn e power modulation can be 
achieved with FACTS devices such as TCSC and STATCOM Excitation s\stem and PSS 
also interact with the torsional modes Combined with control system theorv, eigenvalue 
studies can be used to design controllers of SSR countermeasures Digital time domain 
simulation is most useful tool for the study of transient torque problem In view of the 
above, the objectives and scope of thesis are 

1 To study the subsynchronous transient oscillation damping with the STATCOM 
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2 To propose a new digital control scheme for TCSC control to damp subsynchronous 
resonance oscillations 

3 To study the control interactions between TCSC and povier system stabilizer 

1 5 Outline of the Thesis 

In Chapter 2 the linearized model of generator system with detailed representation of 
stater rotor and mechanical system along with AC network is presented The study 
system is IEEE First bench mark (IEEE FBM) model Eigemalue anahsis is carried out 
with the linearized overall svstem model The results obtained through eigenvalue analysis 
IS validated with time domain simulation and Fourier emalysis The effectiveness of NGH 
damping scheme m mitigating the subsynchronous torque oscillations is studied through 
the detailed time domain simulation using PSCAD/EMTDC 

In Chapter 3 the linearized model of a 12 pulse STA.TCOM along t^ith voltage and 
reactive current controller is presented The control parameters are obtained through 
eigenvalue anah sis of the linearized model The studies are earned out on IEEE FBM 
The results obtained through eigenvalue analysis are validated vuth time domain simulation 
for the study of transient torque oscillations Comparatne stud\ is carried out with NGH 
damping scheme m mitigating subsynchronous torque oscillations 

SSR analysis vith TCSC in continuous time linearized domain is presented in Chapter 
4 Constant angle control of TCSC is investigated The constant angle control is based on 
the philosophy of maintaining the voltage drop across the compensated transmission line 
constant Eigemnlue analysis is carried out to predict the stability of the system and the 
results are validated with time domain simulation 

In Chapter 5, a novel digital control scheme is proposed for TCSC control For this 
purpose a discrete time domain TCSC model is developed The controller is designed using 
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state feedback approach The discrete time model along with the controller is suitably 
interfaced with the generator system model to obtain a homogeneous state space equation 
The stability of the system is anahzed through eigenvalue analysis and the results are 
validated with time domain simulation 

The effect of excitation system that includes power system stabilizer (PSS) on tor 
sional modes, in the presence of TCSC controller is examined in Chapter 6 The control 
interaction study is carried out through eigemalue analysis 

Chapter 7 outlines the conclusions drawn from the thesis and suggests some future 
scope of work 
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Chapter 2 


SSR ANALYSIS OF A SERIES 
COMPENSATED POWER 
SYSTEM 


The stud\ of subsynchronous resonance (SSR) phenomena requires detailed time domain 
simulation of the system considering all relevant nonlmeanties Howeier, the first step 
towards assessing the possibility of occurance of subsynchronous resonance phenomena is 
generalh through the studv of the s\stem in a linearized domain The linearized domain 
analysis is supplemented by time domain analysis for detailed in\ estigation This approach 
is also comenient while designing control for any controllable device m pover svstem 
This hot\e\er requires development of a system model m the Imearized domain For this 
purpose the general practice is to adopt a modular approach in formulating the system 
model with each component being modeled separately to any desired degree of detail The 
various mdividual component models are then suitably interfaced to derive the complete 
system model The development of power system model, following the above mentioned 
approach is described in this chapter Study of subsynchronous resonance phenomena m 
IEEE first benchmark system is carried out using this linearized system model The results 
are validated with detailed t me domain simulation using EMTDC/PSCAD software Also, 
the efficacv of NGH dampi ig scheme [14] in controlling SSR has been studied through 
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digital simulation 


2 1 Power System Component Models 

A single machine infinite bus (SMIB) power system includes generator system and a trans 
mission system The representation of generator system comprising the synchronous ma 
chine and the associated mechanical system is described here along with the transmission 
network These component models are developed in the linearized domain 

2 11 Synchronous Machine Model 

The synchronous machine has three armature windings (a b c) on the stator a field 
winding (/), a damper winding (h) along the direct axis and two damper windings {g and 
k) along the quadrature axis Fig 2 1 shows the layout of the wmdings of the s\ nchronous 
machine In the following analysis, saturation is neglected The model consists of equations 
for stator and rotor windings 
Stator Equations 

The stator is represented bv the two axis equivalent of the three phase windings, a d-axis 
winding on the direct axis and a q axis winding on the quadrature axis The flux linkages 
associated with the d and q axis windings are given bv [4] 

- (1 + Sg) ~ Ratd = Vgd 

7 ( 2 1 ) 

+ (1 + Sg) - Raiq = Vgq 

where ujb is the base rotor speed in rad/sec, Sg is the slip, Ra is the armature resistance, 
Zd and Zg are respectively the d and q axis components of machme current, and Vgd and Vgq 
represent the d and q axis components of the machine terminal voltages T' denotes the 
rate of change of the flux linkage In eq (2 1), all the quantities are expressed m per 
unit 
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To have a common axis of reference with the transmission network, the voltages Vgd 
and Vgq are transformed to synchronously rotating reference frame (D Q reference frame) 
using the following transformation [4] 


'^gd 


cos6 

—s%n8 


VgD 

. '^9Q . 


smS 

cos5 


. ^9Q . 


(2 2 ) 

where VgD and VgQ are the components of voltage along D and Q axis <5 is the angle by 
which d axis leads D axis 
Rotor Equations 


The rotor which has the field winding (/) and a damper winding [g) is represented as [4] 

1 ° (2 3 ) 


do 


where E' and E'^ are d and q axis rotor voltages, and are the open circuit transient 

time constants and Xg are the s\nchronous reactances, X'^ and JT' are the transient 
reactances and Epo is the field voltage referred to the stator side 
The currents %d and z, are related to the flux linkages and voltage sources by 


^d = Vdid + E' 


= Y'Zg - E'^ 


(2 4) 


Substituting Id and Zg from eq (2 4) into eq (2 1) and eliminating Vgd and Vgg using 
eq (2 2), the following equations are obtained after hneanzing the resulting equations 


ATd - - upA^g + ^AE'g - cos^lfgAsg 


+ajBVggA5 - (jJscosSoAvgD + uJssinSoAvgQ 




(2 5) 


:g -usA^d- ^i^A^g- ^^AE'a + UJB^dASg 

-LUBVgdAS - UJBSinSoAVgD - UJBCOSSoAVgQ 

whe e, the subscript 0’ indicates nominal (operating point) values Similarly combining 
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eqs (2 3) and (2 4) and linearizing we get 


AEL = — 

d rpi 

^ go 


^ do L 


X 


1 


Xd 


^-11 A'ir,-^AE,4-A£;^.p 


(2 6 ) 


Defining a state \ector as Ax^i = [A^f^ A'^q AE'^ write eqs (2 5) and 

(2 6) in the follovvnng state space form 


Axsi — [Asi] Axji + [Bsii] Ausii + [Bs 12] AU512 + [Bjia] Au^ (2 7 ) 

where the input \ectors are Aujii = [A(^ AsJ^ Au^p = A^'fz) Aug = [Au^i: Aii^g]^ 
The output equation for stator and rotor variables can be written as 


Ay,i = [Csi] Axsi 


(2 8 ) 


where [C^i] is an identity matrix such that Ay^i = Ax^i The matrices [A^i], [B^n] 
[Bsp] and [B^ia] are defined in Appendix A 

2 12 Mechanical System Representation 


The mechanical s^stem comprises generator, shaft and turbines The shaft can be vieved 
as a mass spring damper system as shown in Fig 2 2 Each major rotating element of the 
mechanical system is modeled as a mass represented b^ its inertia The shaft by a spring 
with its stiffness represented by its spring constant The damping of each mass and shaft 
section is represented by dash pot damping The equation for mass connected by elastic 
shaft sections to mass {i — 1) and mass (% + 1) is given by [4] 


M ^ 4 - D ^ 


+ Arff" + A 1-1 


dt 


di 


+ 


'+1 (f “ ^) + “ '*'+>) = 


(2 9) 




where 

Mj IS the moment of inertia of mass 
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(5i IS the angle of rotation of mass with respect to synchronously rotating reference frame 
A IS the self damping co efficient of mass 

Di i_i IS the mutual damping coefficient between and i — masses 

IS the spring constant of the shaft connecting and i — masses 
Tmi IS the mechanical torque applied on mass 
Te IS the electromagnetic torque on mass 
The inertia Mj is given by 

M, = — (2 10) 

ujb 

where i?, is the inertia constant of mass A.lso, 

5, = ojbSx (2 11 ) 

where s, is the shp of mass For a N mass system, there are N such second order 


D 


12 
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Figure 2 2 Spring mass model of a turbine generator shaft 


differential equations of type of eq (2 9), which an be stacked m the form 


[M] 5 + [D] 5 + [K] <i = [T„] - [Te] 


(2 12 ) 
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where, [M] is the matrix of mass inertia, [D] is the matrix of damping coefficients and [K] is 
the spring constant matrix and Te are the vectors of mechanical and electromagnetic 
torques respectively and 5 is vector of angles of rotation Note that Tg has only one non 
zero element corresponding to the generator rotor Combining eq (2 11) with eq (2 12) 
and then linearizing we get the following state space equation for the mechanical system 
associated with a synchronous machine 

Ax,„ = [A„] Axm + [B,n] Au^ (2 13 ) 


Ay^n — [Ctti AUrn 


(2 14 ) 


where the state vector represents the perturbed values of rotor slips (of different rotors), 
shaft torques and generator rotor angle The input vector u^ represents the perturbed 
value of generator electromagnetic torque applied at the generator rotor mass The output 
vector Yrn consists of perturbed values of generator rotor angle and slip In deriving the 
above equations it is assumed that the mechanical torque input to the turbine is constant 
To illustrate the illustration of derivation of mechanical system equations, let us consider 
the electrical equivalent circuit of a two mass mechanical system shown in Fig 2 3 Mass 1 
is assumed to be a turbine and mass 2 to be a generator In this circuit mass is analogous to 
capacitance, torque is analogous to current and the spring constant is analogous to inverse 
of inductance Following the above procedure the linearized equations for the two mass 
system can be written as 


Asi = ATi — ATi2 — {Di 4* -Dm) ^si + Ass 
2iF2 As2 = — Ar2 + ATi2 — {D 2 + D 12 ) As 2 + X^i2 Asi 


( 215 ) 


Also, the shaft torque (1 e current flowing through the inductance connected between node 
1 and node 2) can be expressed as 


AT12 = Ki 2 (Asi — AS2) 
AJ = (jJq As2 


( 216 ) 
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Figure 2 3 Electrical equivalent of two mass model 

where si is the slip of the mass 1 S2 is the generator rotor slip Tj. is the mechanical torque 
input to the mass 1 and T2 is the generator electrical torque For the two mass system 
shown m Fig 2 3 the state space model is given by 

ZXXjji = [-A-jn] + [Bfn] ^Urrt 

where, Axp, = [ Asi As2 AT12 AJ]^, Ayjn = [A 5 As2]^ and = AT2 

The turbine generator shaft of generating system can be represented as a six mass 
system [ 36 ] It is possible to adopt the above modeling approach to represent the six mass 
mechanical system which consists of a high pressure tuxbine(HP), an intermediate pressure 
turbme(IP), two low pressure turbmes(LPA LPB), a generator(GEN) and an exciter(EXC) 

For the six mass mechanical system the state vector can be written as Ax^ = [Ax„,i Axm2]^ 
where 

Ax,7ii [Asg Asg Asj^pg As^p^ As^p As^p] o,Tid 
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Ax^ 2 — [Ai;, 

\^EXC ^'^LPB^GBN ^'^LPA-’LPB ^Tjp 

-LPA 

The input vector Au^ = ATg and the output vector Ay^ = [A5 A 5 g]^ The expression 
for electrical torque ATg is given b) 


Tg = - '^q'ld 


(2 17) 


Substituting for td ig from eq (2 4) and linearizing eq (2 17) we get the following 
expression 


ATg = 


'vq + E'^ 


x' 


X. 


ATd + 


Tl^d + E'g ^ ^ 


Xa 


A , % 


A^, + + 


X' 


x‘ 


7 


This expression can be rewritten as 


ATg = [C^ei] AXsi 


(2 18) 


where [C„.il = [(*^ - i) 

Axgi = [A^d AT<, AE'd AE'g] 


■ hd±£i + Ik' 


ti 'h. 

^d. 


and as defined earlier 


2 13 Combined Synchronous Machine - Mechanical System Model 

Defining state vector as Ax^ = [Axfj^ Ax^] eqs (2 7) and (2 13) can be combined as 


AX(5 = 


■ Agi 

■A.m 

Axg + 

Bsii 

0 

Augii + 

Bs13 

0 

Aug + 

0 

. '^rn . 


Au„ (2 19) 


From eqs (2 IS) and (2 14), the following is evident Au,n = ATg = CmeiAxgi Auju 
= Ajm = CmA'irn, Using the above expressions, eq 2 19 can be rewritten as 


Axf? = 


Aji BsiiUm 

BAfU'mel 


Axo + 


Bsis 

0 


Aug 


Axg = [ Ag ] Axg + [ Bg ] Aug 


(2 20 ) 
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The output of the generating system model is the currents of the synchronous machine 
defined as Ay^ = [Ai^j Aiq A%n Azq] The currents id and iq are obtained from the 
relationship [4] 

Id _ cosS sin5 td 

iQ \ —sinS cos5 Iq 

Substituting for id and iq from eq (2 4) in eq (2 21) and linearizing we get 

= [Ceil Axg (2 22) 

The matrix [Cd] is defined in Appendix A Combining eqs (2 20) and (2 22), we get the 
output equation for the generating system as 

Aye = [Cfl] Axg + [Dq] Aug (2 23) 

The matrices [Cg] and [Dg] are defined in Appendux A 

2 14 Transmission Network 

The transmission network comprises transmission lines, transformers and capacitors It is a 
general practice to neglect transmission line dynamics in low frequency (0 1-2 Oilz:) electro 
mechanical oscillation studies In contrast the eigenvalue anah sis in the subsynchronous 
frequency range (10 — 402^2;) requires dynamic representation of the transmission network 
based on linearized differential equations In developing the network model, it is assumed 
that the three phases are uncoupled This makes the voltage drops m each phase a function 
of only the current in that phase This assumption simplifies the transmission! network 
equations in the a b c reference frame which are transformed to synchronously rotating 
reference frame {D Q reference frame) The relationship between D — Q reference frame 
and d — q reference frame is given m eq (2 2) A typical series compensated AC network 
IS shown in Fig 2 4 This can be represented m D — Q reference frame as 

Ax// = [A//] Axjv + [B//] Auw 

^Yn — [Cat] Axjv 


(2 21 ) 
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(2 24) 
(2 25) 



where Ax// = [Avcd Avcgf Au^ = [Aigd AzoQf Ay^ = [Avcd Avcof 


An = 


0 

ub 


-OJb 

0 


and Bw = 


* ? 
0 i 


It can be noticed that the input to the transmission network model are the current injections 
obtained as the output of the generator system model Also, the output of the transmission 
network model are the voltages which are the input to the generator svstem model 



Figure 2 4 AC network model 


2 2 Combined Generator-Network System Model 

The state and output equations of the total s\stem can be expressed in the following form 

Ax = [A]Ax+[B]Au (2 26) 

Ay = [C] Ax + [D] Au (2 27) 

where Ax^ = [Ax^ Ax^] Au^ = [^'^g = [Ay^ Ay^] As tl e input 

of one subsystem is related to the output of the other subsystem, the oi ei all input vector 
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can be related to the overall system output vector as 


Au = [F] Ay (2 28) 

The elements of the [F] matrix can be denied from the knowledge of the interconnection 
amongst the various subsystems Now utilizing eq (2 28), eqs (2 26) and (2 27) can be 
combined as 


where 


Ax = [A + BF(J - DF)-^C] Ax = [At] Ax 





[C] 



and [ D ] = 



(2 29) 


2 3 Eigenvalue Analysis of IEEE - FBM 

The system modeling approach described m the previous sections is utilized to study the 
SSR behavior of the IEEE First Benchmark Model [lEEE-FBM] based on eigenvalue analy- 
sis The system details of the IEEE FBM [36] is given in Appendix B The initial conditions 
calculation is given in Appendix C It is assumed that the series capacitor compensates 
70% of the transmission line reactance, i e, Xc = 0 7Xt Mechanical damping is neglected 
m this study The generator delivers 0 9 per unit power at a power factor of 0 9 lag Table 
2 1 gives the eigenvalues of the study system Prom this it can be observed that supersyn- 
chronous and subsv nchronous network mode eigenvalues are stable There are four unstable 
torsional modes The low frequency mode or mode 0 is damped The resonant frequency of 
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Table 2 1 System Eigenvalues 


Series Capacitance = 0 35 
per unit, Xc — 07 Xl 

comments 

4 419 ± j612 43 

3 434 ± jl41 175 

0 0000003 ±j298 176 

0 0019 ± j202 81 

0 0461 ± jl60 27 

0 0373 ± jl27 164 

0 055 ± j99 704 

0 49 ± jlO 56 

3 922 

0 19 

Supersyn Net Mode 
Subsyn Net Mode 
Mode 5 

Mode 4 

Mode 3 

Mode 2 

Mode 1 

Mode 0 


the senes compensated transmission system is obtained from the relationship f^r = 
where L and C are transmission line inductance and capacitance respectively With 70% 
compensation the transmission line resonant frequency is 266 57 rad/s(42 42Hz) and the 
complimentary frequencies are 643 564 rad/s (102 i2Hz) and 110 418 rad/s (17 58Hz) 
From the Table 2 1, it is evident that two of the unstable eigenvalues have frequency close 
to the resonant frequency of 110 418 rad/s (17 58Hz) one at 99 704 rad/s (15 868Hz) and 
another at 127 164 rad/s (20 238772:) This situation may lead to resonance with turbine 
generator shaft To further validate this obsenation, obtained through eigenvalue analysis 
a detailed digital simulation of IEEE FBM is carried out 


2 4 Simulation Studies 


The time domain simulation of IEEE FBM is performed using PSCAD/EMTDG software 
package (version 3 0) to study the response of various shaft torques to a small disturbance 
All system nonlmeanties are represented The field excitation is assumed to be constant 
The s} stem operating cc ndition is same as that taken for eigenvalue analysis A disturbance 
of 0 1 per unit torque is applied at the HP turbine input for a period of 0 1 sec The 
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Time in seconds 


Figure 2 6 Transient torque response for a 3 phase to ground fault with fault impedance 
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(c) Founer spectra of LPA LPB Torque 



Frequency (Tiz) 

Figure 2 7 Founer spectra of LPA LPB torque 
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disturbance is applied at a time t = 0 5 sec Fig 2 o(a) stows the LPA LPB torque as a 
result of this disturbance The LPA LPB torque in the mterval 4 0 sec to 4 1 sec is shown 
in Fig 2 5 (b) The Fourier spectra of the LPA LPB shaft torque shown in Fig 2 5 (c) 
This indicates a predominant presence of 20 &2Hz component which is close to the mode 
2 frequency obtained in eigenvalue analysis 

To study the system behavior under large disturbance the study system is subjected 
to a three phase to ground fault through a fault impedance of 0 04 per unit at infinite bus 
end The fault is applied at 0 5 secs and is cleared after 4 5 c\ cles The shaft torques are 
shown Fig 2 6 

From the plots given Fig 2 6 it can be observed that the shaft torques increase with 
time as a result of torque amplification This behavior is mdicative of the unstable system 
operation due to the onset of SSR phenomenon A s imil ar observation regarding unstable 
system behavior is evident from the eigenvalue results gi\en m Table 2 1 The LPA-LPB 
torque of Fig 2 6 is reproduced m Fig 2 7(a) Its -variation over the interval 4 0 sec to 
4 1 sec IS sho-wn m Fig 2 7 (b) The Fourier spectra of LPA LPB torque is shown m Fig 
2 7 (c) This plot also shows the presence of 20-ffz component, a frequency which is close 
to mode 2 frequency obtained in eigenvalue analysis 

From the time domain simulation, it is evident that the amplitude of shaft torques 
increase which mav result in the turbine generator shaft damage Many generator based 
counter measures are suggested m the literature to deal with transient torque problem 
[4, 11, 14,37,38] One of the capacitor based counter countermeasure to suppress transient 
oscillations due to subsynchronous resonance is NGH dampmg scheme [14, 15] In the 
following section the application of NGH dampmg scheme to IEEE first benchmark model 
is discussed 
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2 5 Study of lEEE-FBM with NGH Damping Scheme 

The basic NGH scheme for one phase includes a linear resistor in series with anti parallel 
thyristors connected across the capacitor an shown in Fig 2 8 The anti parallel thyristor 


SERIES 

CAPACITOR 



CAPACITOR 

VOLTAGE 

SIGNAL 


THYRISTOR 

CONTROL 


Figure 2 8 Basic NGH scheme 


pair in the NGH damper functions as a bidirectional switch The gating of the pair is 
so controlled that the current can flow in either direction depending on the direction of 
the line current For normal operating condition, successive zero crossing points of the 
capacitor voltage are timed In the steady state this time will be equal to half cycle period 
of the system frequency (j, where T is the time period of the voltage) If and when half 
cycle period deviates from the set time of the thyristor switches are closed to discharge 
the capacitor energy through the linear resistor 

The system under study is IEEE FBM with a series capacitor of (0 35 per unit) reac 
tance which compensates 70 % of the transmission line reactance The series capacito is 
assumed to be made up of two banks, one having 0 301 per unit capacitance and the otiier 
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0 049 per unit capacitance The NGH scheme is applied to 0 049 per unit capacitor bank 
In the study all nonlineanties are represented and the field voltage is assumed constant 
Self damping of 0 2 per unit and mutual damping of 0 3 per unit is added to the haft To 
study the effectiveness of NGH damping scheme in mitigating SSR, the following studies 
are carried out 

(i) Capacitor insertion test 

(ii) Three phase short circuit test 
(ill) High speed reclosure test 

(i) Capacitor insertion 

During capacitor insertion operations, subsvnchronous frequency oscillations are present 
due to the interaction between the line inductance and series capacitors and these oscilla- 
tions mav excite the torsional modes A capacitor insertion study is carried out without 
NGH device In this study, the generator is delivermg 0 9 per unit power Imtially the 
capacitor segment of 0 301 per unit is in semce and the other capacitor segment of 0 049 
per unit is inserted at time t = 0 6 secs Fig 2 9 shovs the transient torque response with 
the fixed capacitor insertion It can be seen that durmg the fixed capacitor insenion the 
GEN E\G torque is excited and the oscillations are growing as a result of subsynchronous 
resonance phenomenon The NGH device is expected to suppress these oscillations Con 
sequentlv the capacitor segment of 0 049 per unit is inserted with NGH device at t = 0 5 
sec, while the fixed capacitor segment of 0 301 per unit is considered to be in semce The 
GEN E\C torque and LPA LPB torque responses axe shown m Fig 2 10 From the Fig 
2 10 it IS evident that the NGH device is effective m stabilizing the shaft torques durmg 
capacitor insertion by damping the subsynchronous oscillations 

(ii) Three-Phase short circuit fault 

In this study the generator is deliveimg 0 9 per unit power Three phase to ground faults. 
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Figure 2 9 Transient torque response with fixed capacitor insertion 

with, duration of one cycle (0 01666 sec) and 4 5 cycles (0 075 sec) are applied at the infinite 
bus end through a fault impedance of 0 04 per unit The fault is applied at time t = 0 5 sec 
The shaft torques response IS shown in Fig 2 11 and Fig 2 12 From the responses shown 
m Figs 2 6, 2 11 and 2 12, it can be observed that with the NGH damper the transient 
oscillations do not grow and on the contrary have gradually decreasing amplitude This is 
indicative of stable system operation The Gen Exc torque oscillations are reduced faster 
compared to the other shaft torques The peak values of the different shaft torques for the 
two different fault durations are given in Table 2 2 

Table 2 2 Peak shaft torques for different fault clearing times 


Fault clearing 

] 

peak value of torques in per unit 

time 

Gen Exc 

LPB Gen 

LPA-LPB 

IP-LPA 

HP IP 

m seconds 

torque 

torque 

torque 

torque 

torque 

0 016 

0 456 

2 751 

2 846 

1 838 

1 054 

0 075 

0 4248 

2 465 

2 50 

1 52 

0 819 
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LPA LPB Torque in per unit 



Figure 2 10 Transient torque response with insertion of capacitor having NGH damping scheme 
( 111) High speed reclosure 

When a generator is connected to multiple lines and a fault occurs close to the generator 
on one line, the system voltage drops and the generator is no longer able to deliver the 
prefault power to the network Consequently the rotor angle increases The restoration 
of the voltage at the instant of line tripping at increased load angle causes stress on the 
shaft segments The reclosure of the line has most significant effect on the shaft transient 
torques A successful reclosure results in high transient torques due to the sudden return 
of load An unsuccessful reclosure subjects the system to another fault In the case of a 
single tied generator, the load angle increases after fault occurance, during short circuit 
and the reclosing time intervals A successful reclosure with a high magnitude of phase 
difference between the generator and the system causes severe stress on shaft segments 
Thus the reclosing action for a single tied generator results in higher shaft torques in series 
compensated systems In the present study, the generator is delivering 0 5 per unit power 
to the infinite bus A three phase to ground fault is applied at t = 0 5 sec The fault is 
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cleared after 3 cycles by the breakers at both ends of the line The breakers are closed 
successfully in 15 cycles after fault clearance Fig 2 13 shows the dynamic response of 
shaft torques to a 3 15 cycle successful reclosure test and it can be observed that, even 
at a reduced power level high speed reclosure test is producing peak torques greater than 
those encountered during the three phase fault case shown in Fig 2 11 The NGH scheme 
IS able to render the damping effect on shaft torque oscillations 

2 6 Conclusions 

In this chapter the torsional mode stability information obtained through eigenvalue anal 
ysis using linearized system model is verified with detailed time domain simulation The 
control of transient torques with NGH dampmg scheme is illustrated with time domain 
simulation The simulation results pertainmg to capacitor insertion case reveal that the 
NGH damping scheme is able to stabilize the transient torque oscillations, but the dampmg 
effect is inadequate For the case of three phase fault disturbance the damping effect is 
found to be more pronounced on Gen Exc torque compared to the other shaft torques 
The rate of damping of oscillations of LPA LPB and LPB Gen torques is moderate This 
mav be due to the chosen value of capacitance Though NGH dampmg scheme arrests 
the growth of transient oscillations but its effect on peak torque reduction is not signifi 
cant The NGH device do not have the capability of increasing the power transfer and its 
function IS limited to dampmg the transient torques 

Modulation of real power and reactive power are the effective means of damping the 
power system oscillations and this modulation can be obtained by incorporating FACTS 
devices such as TCSC and STATCOM In the next chapter, the control of subsynchronous 
resonance phenomena with the STATCOM is discussed 
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Figure 2 11 Transient torque response for a 3 phase fault (4 5 cycles) with NGH scheme 
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Chapter 3 

SSR ANALYSIS WITH STATCOM 


The Flexible AC Transmission System devices will need to be used along with the comen 
tional series capacitors to achieve the objective of impro\ements in the loadability of the 
system and extended performance The availability of high power gate turnoff thyristors 
(GTO) have recently led to the development of a new equipment called Static Compen 
sator (STATCOM), which is a second generation FACTS device based on voltage source 
inverters STATCOM uses self commutating devices like CTOs and is an advanced form 
of SVC The operating and functional characteristics of STkTCOM are, however different 
from those of S\ C The STATCOM does not use capacitor or reactor banks to produce 
reactive power It operates as a controllable electromc synchronous voltage source con 
nected to the line through the leakage inductance of a coupling transformer and draws 
reactive current from the line in essentially the same wa\ as a rotating Synchronous Con 
denser (without the mechanical inertia and long excitation time constant) The controlled 
synchronous voltage is maintained in phase wnth the line voltage It can, thus, draw either 
inductive or capacitive current by varying the magnitude of this voltage Furthermore, the 
DC capacitor bank, to which the voltage source inverter is connected, does not play an 
active role m the VAR generation 

The subsynchronous resonance oscillations can be damped by modulating the real and 
reactive power The reactive power modulation based countermeasure is considered in this 
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chapter, where we consider the STATCOM as a reactive power source The STATCOM is 
connected at the generator terminal and operated m voltage control mode 

3 1 System Modeling 

The s\stem under study is the lEEE-FBM with fixed series compensation The config 
uration of the study system which is shown in Fig 3 1 consists of a generator feeding 
power to an infinite bus through a series compensated transmission line The STATCOM 
IS connected at the generator terminals The system can be divided into the following 
subsystems 

(i) Generator System 

(ii) AC Network 
(ill) STATCOM 

The generator system model and the network model are the same as given in Chapter 
2 For the sake of completeness the models are reproduced below 
Generator System 

The state space model of the generator system is gi\en by 

Axq = [Ag] Axg + [Bg] Aug (3 1) 

where Ax^ = [Axfj^ are the state vectors of stator and mechanical 

system, Aug = [^Vgo Au^q]^ The output equation is given by 

Ayg = [Cg] Axg + [Dg] Aug (3 2) 

where Ayg = [Azg A?q Aig A iq] The matrices [Ag], [Bg], [Cg] and [Dg] are defined 
m Chapter 2 
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AC Network 


The state and output equations of the AC network are given by 

Ax/v = Ax;v + Pw] Au;\^ (3 3) 

Ayat = [Cat] Ay.^ (3 4) 

where state vector contains D Q axes capacitor voltages i e Ax/y = [Avcd AvcqY' 
The input vector contains D Q axes line currents and D Q axes STATCOM currents i e 
Aua^ = [Azgd Aiqq Aisd Aisq]^ The output vector is given as Ajn = [Avcd Avcof 

3 11 STATCOM Representation 

The STATCOM is based on the voltage source inverter (VSI) using Gate Turn Off (GTO) 
thyristor valves A basic 6 pulse VSI circuit is shown in Fig 3 2 This consists of sl\ GTO 
thyristors each having an anti parallel diode and a capacitor connected on the DC side 
The VSI IS connected to the system bus through a small reactance which is the leakage 
reactance of the coupling transformer 

The magnitude of the VSI output voltage can be controlled to produce a fundamental 
output voltage m phase with the A.C system voltage If the VSI output -voltage is greater 
than the AC system voltage, capacitive current is drawn from the power system and con 
versely, if the VSI output voltage is lower than the AC system voltage, inductive current 
IS drawn By varying the phase angle of the inverter output voltage relative to the AC 
system voltage the output voltage of the STATCOM is varied thereby varying the reactive 
power supplied or drawn Also, the STATCOM can negotiate the exchange of real power 
with the AC system by varying the phase angle of the inverter output voltage That is, the 
STATCOM supplies real power to the AC system from its DC capacitor if the STATCOM 
output vol age is made to lead the AC system voltage By the same token, the STATCOM 
absorbs re il power from the AC system, if the STATCOM output voltage is made to lag 
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Figure 3 2 A basic 6 pulse voltage source inverter 


the AC system voltage This particular property is used to maintain the capacitor voltage 
at the desired level 


3 12 Modeling Of STATCOM 


Consider the circuit of Fig 3 2 A set of three phase voltage is produced from the input 
DC voltage The GTOs are assumed to have negligible turn off and turn on time For the 
discussion presented below, we shall neglect both AC and DC side harmonics and consider 
only the fundamental frequency components The dynamics of a six pulse STATCOM can 
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where v’ vl and ti* are the three phase instantaneous voltages at the inverter AC side 
terminals R and are the resistance and leakage reactance of the transformer, Rp rep 
resents the switching losses m the in\erter Vdc and idc are the voltage across and current 
through the capacitor and Zc are the instantaneous currents on the AC side and Va 
Vb and Vc are the instantaneous voltages of the AC bus The instantaneous three phase \C 
bus phase voltages are given by 

Va — VmSin {ulst + 9) 

Vb = VmSin (uJBt + 9 — 120°) (3 6) 

■Wc = VmSin (ust + 9 + 120°) 

where where 9 is the AC bus voltage angle defined with respect to the infinite bus voltage 
of Pig 3 1 and 

For a SIX pulse inverter with pulse duration of 60° the peak value of the fundamen 
tal component of the phase voltage is ^Vdc consequently, the three phase instantaneous 
voltages at the inverter AC terminals can be expressed as 

2 

Va = Vdc-szn {ust + /3 + 9) 

TT 

vl = Vdc— sm (uijBt + j3 + 9 — 120°) (3 7 ) 

TT 

2 

V = Vdc— sin + /d + 0 + 120°) 

TT 

where /? is the angle between the imerter output voltage and AC bus voltage and ujb is 
the base frequency 

Assuming power balance between the AC and DC sides, the following equation can be 
written inverter [31] 

Vdc^dc - v'ala + vllb + (3 8 ) 
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Defining 7 = {wat + /? + 0) we write the dc current expression using eq (3 8) 


2 

idc = - [^QS^n (7) + Zfrszn (7 ~ 120°) + icSin ( / + 120°)] (3 9) 

Substituting eqs (3 6) and (3 7) in eq (3 5) we get 
u)b **2 

= --TT VdcSin (7) - UmSin (wat + 9) 

As I TT 

lb = — — Rib — —VdcSin (7 — 120°) — VmSin {ust + 9 — 120°) (3 10) 

-^S ^ ^ 

cj ' 2 

2c = VdcSin (7 + 120°) — VmSin (ust + 9 + 120°) 

As L TT 

The above equation can be written in the following state space form 

Xabc = [A] Xabc + [Bi] V^c + [Ba] Ua6c (3 H) 

where = [la ib ic], = [va Vb Vc] 

Substituting eq (3 9) in eq (3 5), the DC capacitor voltage cau be written as 

Vdc = -^ ^ + -(2a52n(7) + 2i,s2n(7 - 120°) +2cS2n(7 + 120°)) (3 12) 

be iRp TT 

The quantities in the abc reference frame can be expressed in DQO reference frame using 
the following transformation 

■ 2a 1 r cos(ajat) siniuBt) 75 1 [ 

cos(u;at - 120°) S2n(u;Bt - 120°) ^ iq = [T] [2aQo] 

_ 2c _ COs{uJBt + 120°) Sm{(jJBt + 120°) . 

(3 13) 

Premultiplying both sides of eq (3 11) by T~^, the following expression is obtained 

[T-'jXatc = [T-'] [ 4]Xa6c + [T"'] [5i]vdc + [T"'] [52]Ua6c (3 14) 

E quation(3 14) can be rewritten as 

[T-']xa6c = [T-i][A][T][T-^]xa6c + [T-'][5i]vdc + [T-^][52]Ua^ (3 15) 
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Using eq (3 13) the following expression is written 


T ^Xaf,c = ^DQO + T ^Xaic = ^DQO + T ^TXdqq 


1-1 


(3 16) 


Substituting eq (3 16) in eq (3 15) gives 


xnqo = [^] XjPqo 


T“1T 


^DQO + [-^l] Vdc + T ^ [^2] MdQO (3 17) 


where x^qq = [t^D Isq ^so], '^dqo = Kr? '^sq ^so] the quantities Isd, T'sQ, ho and Vsd Vsq Vso 
are the STATCOM currents and voltages in DQO reference frame 

Prom eq (3 17) the DQ components of STATCOM currents can be obtained as 

Rlob sin((3 + 9) ojb 

'>-sD — — ^B't'sQ + ^B'^'^dc ^7 ~'^sD 


X. 


A, 


X 


Rjuob , , cos{P + 9) ujb 

IsQ = OJbIsD + n^s/CUtic W 

A e A fl 


‘•S Xj 


Similarly, using the quantities expressed in DQO reference frame eq (3 12) can be written 
as 

, s%n{P + 9) cos{P + 9) ojb 

Vdc — '^COsk r IsD ^sk r ^sQ - _ 

Oq OfH; 

where A: = ^ for six pulse converter and k = ^ for a 12 pulse converter 
The above expressions for IsDi'I'sQ and Vic can be combined as 


X5S = [Ao] 'yisE + [Bo] Us 


(3 18) 


where xsb = [isi? %sQ Vdcf', us = and [Aq] = [J] The elements of the [J] matrix 

IS given below 

T R^b j t j^sin{l3 + 9) 

= 77“) >^12 = ~n;j5, ^13 = UIbK- 


X 


X 


T T _ r _ ^cos{l3 + 9) 

J21 — U)B, J27 — 7r ~) — W ^/ C - 


X 


X 


sin{^ + 9) cos{P + 9) CUB 

J3I = —UbK r J32 = —OJsk Jz 3 = 


i>cRp 
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The matrix [Bq] is given by 


[ Bo ] 


X 


0 

0 


0 

X 

0 


'VsQ 


It may be noted that STATCOM voltage Vs is same as the generator voltage Vg since both 
STATCOM and the generator are connected to the same bus as shown in Fig 3 1 The 
angle 9 can be expressed as 


9 = tan ^ (3 19) 

k'^sqJ 

Linearizing the eq (3 19) about an operating point, we get 

- ^Av,q (3 20) 

where is the STATCOM bus voltage at the nominal operating point Linearizing eq 
(3 18) and utilize eq (3 20) the following state space equation is obtained 


Axse = [Ass] Axss + [Bssi] Aus + [Bs£2] A/3 


(3 21) 


where Ass = Aq, Bssi = G, and Bssz = H The elements of the [G] matrix are given 
below 

k Vdco cos{Po+9o) VsQo' 


G 


11 


'x, 


1 - 




tjJB VsDO , ra , a \ 

! Gi 2 — — — 77 — VdcO COs(Po + ^0)) 


G 21 — -^k Vdco s%n{l3o + 9o), G 22 = 


k VdcQ sm(Po + do) VsDO 


^sO 


Gn = kci ^32 = -kci and Ci = cos{Po + 9o) + %sqo sm(/?o + 9o) 

The elements of the [H] matrix are given below 
Hn =iVBk Vdco H 2 x = ub k Vdco and 

Hzx = (^sjDo cos{^o + 9o) — %sQo sm[pQ + ^ 0 )) 

The output equation is given by 


Ay SB 





0 

0 

1 

l^lsQ 


_ 0 

1 

j 

0 
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Axss = [Ose] Axss 

fEmiAL LfBUABY 

I I 7 KANFMI 


(3 22) 


A 13X0L66 
















The STATCOM controller dynamics can be expressed through the following equations 


3^1 = + ^{vs- KdIr) 

^2 ” (VTsfgtat ^l) 

X3 = ^(-Xs+Ir) 

TTT 


(3 23) 


Xii — K 2 {X2 T Kp {^refstat ^l) ^s) 

The reactive power Qg of the STA.TCOM is expressed in terms of D — Q components of 
voltage and current as 

Qs ~ '^sQ'^sD Xg£)lsQ 

From the above equation we can wTite the expression for STATCOM current as 

VsQlsD ~ XsD'l'sQ 


iR 


(3 24) 


In the linearized domain eq (3 25) can be written as 


[VsQqAIsD — XsDQ^'I'sq] + 


+ 


UDO __ 1 
t 0 


T ^^<30 






An 


(3 25) 


sQ 


and the STATCOM bus voltage can be expressed in its D - Q components as 




It may be noted that Vg is same as I g shown m Fig 3 1 Linearizing the above equation we 
get 

(U 5 O + AUs)^ = {VgdO + AVsd)^ + {fJsQO + 

Neglecting the higher order terms, the above equation can be approximated as 

A«. = '^Av.v + (3 26) 


VgO 


n.,0 


Combining eq (3 23) wuth eqs (3 25) and (3 26), we get 


= [-A-sc] Axsc + [B5C1] Ausci + [6502] Au5C2 


(3 27) 
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where Ax|’(^ = [Axi Arcj A^s AX4] = [Az^p Az^q] and Au^^j = [Az^jd Au^q] 

The output equation is given by 


Aysc = A/3 = [Csc] Ax^c (3 28) 

The matrices [A 5 C] [B^ci] [B 5 C 2 ] and [Csc] are defined in Appendix D 

The input to the STA.TCOM controller state space model is the D — Q components 
of STA.TCOM current which is obtained as the output of the STATCOM model The 
output of the STA.TCOM controller is A^ which is the which is one of the input in the 
STATCOM model Based on this mter component relationship, the STATCOM controller 
and the STATCOM models can be combined using eqs (3 21), (3 22), (3 27) and (3 28) as 

Axs = [A 5 ] Axg + [B 5 ] Aug (3 29) 


where. 


[ AAs ] = 


Asfi B5E2 Csc 
BsciCs£ Asc 


and Aus = Ausc 2 Also, 




Bs£1 

Bsc2 


Ays = 


AZ5£) 

Az^q 

Azso 

AZsQ 


= [Cs] Axs + [Ds] Aus 


(3 30) 


The matrices [Cs] and [Ds] are defined m Appendnx D Equations (3 29) and (3 30) give 
the state space model of the STATCOM including its controller 


3 2 Combined Generator-Network-STATCOM Model 


The various subsvstem models given by eqs (3 1) to (3 4), (3 29) and (3 30) can be combined 
to get the overall system model The interconnection diagram is shown in Fig 3 4 Prom 
the diagram it can be seen that the input variable of STATCOM model are related to 
the output variables of the network model The network input \anables are obtained from 
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the generator system and the STATCOM models The various subsystem models can be 
interconnected m the manner 

Au = [F] Ay (3 31) 

where Au = [Auq Aug Auj^-]^, Ay = [Ay^ Ayg Ay|)] and Ax = [Axg Axg Ax^] 
The elements of matrix [F] provide the interfacing between various component models as 
shown in Fig 3 4 


Generator System AC Network 



STATCOM 


Figure 3 4 Interconnection pattern 

For the AC network the injected current is the algebraic sum of generator current 
and STATCOM current Hence the input vector in the network model will change to 
accommodate the STA COM output variables The input and output vectors of individual 

subsystems are gl^ en bt lo\s 
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For the generator 


Aug = [l^VgD /^Vgof and Ay a = [Aid Azq Aid Azq]^ 
For the AC network 


Auiv = [Aiqd Aiqq At,D A Isq]'^ and Ay at = [Avcd 
For STATCOM 

Au5 = [Aus£) Avsof and Ays = Aisq Aisd Ai^q]^ 

The state and output equations of the component models are combined to obtain the 
state and output equation of the entire system m the form 



Ax = [A] Ax + [B] Au 

(3 32) 


Ay = [C] Ax + [D] Au 

(3 33) 

where 

[A] = diag [Ag A5 A^] [B] = diag [Bg Bs Bj,] , 



[C] = diag [Cg C 5 Civ] and [D] = dzag [Dg D 5 Djv] 


Substituting eq (3 33 ) m eq (3 31) to eliminate Ay and incorporating the resulting ex 

pression for Au m eq (3 32) gives 


Ax = [A + BF(I - DF)-^C] Ax = [At] Ax 

(3 34) 


The details of F matrix are given m Appendix D 
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3 3 Eigenvalue Analysis 


To illustrate the development of the system model the eigenvalue analvsis of IEEE FBM 
IS under taken In the IEEE FBM the generator is delivenng 0 9 per unit pover to 
transmission system at 0 9 power factor lag The generator excitation is assumed to be 
constant Self damping of 0 2 per unit and mutual damping of 0 3 per unit is added to 
the shaft to account for steam damping [13] The value of senes capacitor is chosen to 
compensate 70% of the transmission line reactance(i e Xc = 0 3o per unit) ST-VTCOM 
data IS given m Appendix D A. 12 pulse STATCOM is assumed to be supplying reactive 
power of 0 14 per unit The initial condition calculation of STATCOM is given m Appendix 
D The overall system eigenvalues are given in Table 3 1 

The study is initially carried out considering that the STATCOM controller has a 
voltage regulator only The voltage controller gams are taken as iCp = 0 74 and Kj = 10 
To neglect the effect of the reactive current controller the gains Ki and Ki are set to zero 
From the system eigenvalues given in Table 3 1, it can be seen that without the STATCOM 
the mode 2 is unstable Introduction of the STATCOM with onh the voltage results m the 
stabilization of mode 2, while at the same time it destabilizes the mode 4 The damping 
of mode 3 and mode 1 is increased The networks mode damping is also unproved 

In a subsequent study, both the voltage and the reactive current regulators are con 
sidered The reactive current controller gains are chosen as iFi = 0 13 and JCo = 0 01 
From the eigenvalue results it can be observed that STATCOAI increases the dampmg of 
all the torsional modes when both the controllers are active The dampmg of the network 
modes is greatly enhanced Compared to the previous case, the dampmg of mode 2 is 
increased, whereas the dampmg of mode 1 and mode 3 is decreases The subsynchronous 
network mode clamping is marginally decreased, whereas the supersynchronous network 
mode dampmg s increased The frequency of oscillation of subsynchronous network mode 
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IS decreased compared to that in the earlier case The effect of the STATCOM in improv 
ing the damping of torsional modes as observed from the eigenvalue analysis is further 
validated through dynamic digital simulation 


Table 3 1 Eigenvalues for different control strategies of STATCOM 


Without 

STATCOM 

With STATCOM 
voltage control 
alone 

With STATCOM 
voltage control and 
reactive current 
control 

comments 

4 42± j612 43 

10 17±63o 19 

13 09±637 92 

Supsvn Net Mode 

3 45± jl41 11 

28 93± jl69 79 

28 04± jl34 69 

Subsvn Net Mode 

1 38 ± j298 16 

2 05±j298 15 

2 05±j298 15 

Mode 0 

0 28 ± j203 03 

0 65 ± j204 88 

0 534 ± j204 24 

Mode 4 

0 89 ±jl59 12 

1 56 ± jl60 72 

1 27 ± jl60 11 

Mode 3 

0 03±jl29 14 

0 217±jl2919 

0 218±jl29 128 

Mode 2 

0 ll±j98 6ol 

0 158±jl00 159 

0 88±j99 64 

Mode 1 

0 43±jl0 56 

0 71±j26 78 

278 941±4189 75 

250 

250 

2 24±jl8 44 

285 34±4192 48 
156 13±253 33 

0 32±0 35 

250 

Mode 0 

3 91 

-3 695 

35 


0 20 

0 168 

83 47 

0 079 



3 4 Time Domain Simulation 

Time domain simulation is carried out to observe the shaft torque response durmg major 
disturbance using PSCAD/EMTDC package version 3 0 It may be mentioned here that 
the torque amplification is severe at higher compensation levels [36] The chosen capacitive 
compensation provides suitable test case for torque amplification All system nonhneari 
ties are incorporated in detail Dynamics of the excitatio i system is neglected Governor 
action is not considered The STATCOM is considered to be equipped vith both voltage 
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and reactive current controllers The generator is delivering 0 9 per unit power A three 
phase to ground fault of 4 o cycle duration with 0 04 per unit fault impedance is applied at 
the infinite bus at time t = 0 5 sec The plots of the shaft torques are shown in Fig 3 5 It 
can be observed from the results that the STATCOM is effective m damping the torsional 
oscillations For the same level of transmission line compensation the peak shaft torques 
associated with NGH damping scheme and STATCOM are tabulated m Table 3 2 It can 
be seen that there is a significant reduction in the peak values of GEN EXC and LP4. LPB 
torques with ST-kTGOM The reduction in peak values is marginal in the case of IP LPA 
and HP IP torques Though there is marginal increase m peak value of LPB- Gen torque, 
the amplitude of oscillations are progressively decreasing It can be seen from Fig 3 6 that 
the STATCOM is also able to regulate the generator terminal voltage effectivelv It may 
be mentioned that no attempt has been made here to optimize the controller parameters 
of STATCOM 


Table 3 2 Peak shaft torques with different devices 


Name of the 

] 

peak value of torques in per unit 

Dev ice 

Gen Exc 

LPB Gen 

LP\ LPB 

IP LP4. 

HP IP 


torque 

torque 

torque 

torque 

torque 

NGH damper 

0 4248 

2 465 

2o0 

1 52 

0 819 

STATCOM 

0 3085 

2 5 

2 16 

1 47 

08 


3 5 Conclusions 

In this chapter the effectiveness of the reactive power modulation feature of STATCOM 
in damping the subsynchronous torsional oscillations is studied A linearized model of 
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the STATCOM along with the IEEE FBM is derived to study the transient torque 
problem The effect of both the \oltage and reactive current controllers of STA.TCOM 
IS investigated The voltage control scheme alone is not able to provide entire damping 
to torsional oscillations The voltage control along with the reactive current control is 
able to damp the torsional oscillations The control scheme is implemented usmg local 
measurements alone 

The Shunt compensators provide an indirect approach to mcrease the pover transfer 
over a transmission line to certain extent The series capacitors provide direct approach 
to mcrease the power transmission capability which is very often the most economical 
solution In the above discussion the damping of SSR oscillations is demonstrated using 
reactive power modulation The real power modulation can also be used for damping SSR 
oscillations The TCSC is one of the FACTS devices that can be used for real pover mod 
ulation In the next chapter the damping of SSR phenomenon using TCSC is presented 
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Chapter 4 

SSR ANALYSIS WITH TCSC 


In the previous chapter SSR analysis with the STATCOM is discussed The th\ristor 
controlled series capacitor (TCSC) is one of the FACTS device considered for control 
of power flow over the transmission line In this chapter, SSR analysis with TCSC is 
presented A simple proportional mtegral control is used for TCSC control The TCSC 
model IS same as that given in [21] The study is carried out using eigenvalue anal-^sis and 
the results are validated through PSCA.D/EMTDC simulation 

4 1 System Modeling 

The system under studv is the IEEE FBM with a controlled series compensation as shown 
in Fig 4 1 The variable reactance Xc shown in Fig 4 1 comprises of fixed capacitor and 
a TCSC It can be seen that the studv system is the same as that given in Fig 2 4 except 
for variable Xc Hence the eigenvalue analysis technique given in Chapter 2 is directly 
applicable here However, to facilitate the inclusion of the continuous tune TCSC model 
and its associated PI controller, the AC network equations will be different from those 
given in Chapter 2 The turbine generator model is also the same as given m Chapter 2 
For the sake of completeness, the turbine generator model is reproduced here 

Ax(5 = [Ag] Axg + [Bg] Aug (4 1) 
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Figure 4 1 IEEE FBM -with Controlled Senes Compensation 

where Ax^ = [Ax^i Ax^], Aug = [Ati^g Avgof, Axgi and Axm are the state vectors 
associated with the S'vnchronous machine and the mechanical system The output equation 
of the generator system is given by 

Aye = [Cg] Axg + [Dg] Aug (4 2) 

where Ay a = [Az^ Aiq Ain Azq] The matrices [Ag] [Bg] [Cq] and [Dg] are defined 
m Chapter 2 

4 11 AC Network Representation 

The AC network consists of transformers, transmission line a fixed capacitor and a TCSC 
and it can be represented as a two port network with the generator current entering into 
one port and the other port being connected to the infinite bus This is shown in Fig 4 2, 
where, zg represents the current through the line, ej represents the infinite bus voltage 
and vcF represents the voltage across the fixed capacitor Cp It is assumed that machine 
output currents in and iq are equal to the network injection currents, i e , tan and igq 
respectively 

The TCSC model given in [21] is developed based on the representation of voltages 


65 




Figure 4 2 AC Network 

and currents as time varying Fourier series quantities and utilizing the short term Fourier 
coefficients The Fourier senes is truncated to get the fundamental components of voltages 
and currents in a quasi static phasor representation The voltage and current phasors 
are used to describe TCSC dynamics with the assumption that line current is essentially 
sinusoidal The quasi steadv state model of the TCSC is mentioned belov and the details 
are given in Appendix E 

CtCSC = % + y (rr) 

^TCSCe.ff\^) 

where, XTcsceff{<^) is the effective TCSC reactance and cr is the conduction angle The 
conduction angle is related to firing angle by the relation cr = 7r~2a The network dynamic 
equations are 

Cp VCF — I'G (4 3) 

CtCSC Vrasa = + vr^- 7- T 

■^TCSC£ff[G) 


(4 4) 
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The network equations m D — Q domain are given by 


'^CFD -p:, ^CFQ 

Uf 

'^GQ 

‘^CFQ — ^ (jJb ^CFD 

Uf 


'^GD 


^TCSCD 


^ '^tCSCQ /O \r /\ 

UtCSC UtcSC ^TCSC\^) 


^TCSCD 


'^GO 


+ U)f '^tcscd 


''rcscQ 


^TCSCQ ^ ^TCSCD ^ n V* (.r\ 

Ur CSC UtcSC ^TCSC[or) 

These equations are linearized to get the following state space form 


Ax^ = [Aiv, Axjv + [B//_ Auiv 


(4 5) 


where Axjv — [Aucfi? ^‘^cfq ^'^tcscd ^'^tcscq^ Au^r = \_£\iaD Aiqq Act Aiqd A?gq] 
The generator terminal voltage is given by 


— fib + ^CF + '^tgsC 1'G + T'G 


The above equation is transformed mto D Q reference frame as 


VgD = SbD + VCFD + ^tcscd + '>‘GD + Iqd + Wb Ti iQQ 

'^gQ = fibQ + ^CFQ t ^tcscq + ^F 'I'GQ + -^L ^GQ — FJb Li Iqd 

where Ri and Li are the resistance and inductance of the transmission hne respectively 
Linearizing the above equations and combining with eq (4 5) we can i\Tite the output 
equation in the following form 

— [Cw] Axjv + [L)n] Amn (4 6) 

where Ay^ = [Avgo AvgQ Au^csco ^'^tcscq] matrices [A^v], [Biv] [C/v] and [D^] 

are defined in Appendix E 
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4 12 TCSC Controller 


The TCSC is represented as a vanable capacitive reactance in series with the transmis 
Sion line By changing the senes compensation level using automatic controller (affected 
through the use of thyristors) the pover flow through the line can be controlled to improve 
both steady state and transient stability and to damp out the oscillations The purpose 
of the controller is to maintain the TCSC reactance at the desired value specifled by a 
the flrmg angle of thyristors The TCSC reactance remains constant in steadv state and 
is modulated during transient conditions in response to the control signal The objective 
of the control here is to maintain the machine angle constant Assuming that both infinite 
bus and generator voltage being regulated, this amounts to keeping the voltage across the 
line constant Let this voltage across the line be denoted by V , and the current flowing 
through the line by Iq, line mductne reactance by Xl and fixed capacitue reactance by 
XcF we, thus, have the following requirement 

|V| = \Ig{Xl - \cf) - Vtcsc\ = constant = Ki (4 7) 

In the above equation, the hne resistance is neglected and the line current and the TCSC 
voltage are represented bv their rms ralues For the study system the increased power flow 
through the compensated line is associated with the increase in the current Iq through the 
hne In order to maintain the voltage drop V across the line at a constant value the above 
equation implies that Xtcsc should mcrease Thus, m this case, the senes compensation 
must increase with the total power dispatch so that a higher power is transferred at a given 
angle This is called constant angle strategy 

The TCSC control system is shoivn m Fig 4 3 [22] The TCSC voltage is passed 
through a gam Ky = x^-x^ connected to line current reference junction The line 

current %g is monitored and compared to a desired reference current The algebr ic sum of 
the reference current, line current and the signal proportional to TCSC \oltage is an error 
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signal that is fed to proportional integral (PI) controller The output of the PI controller 
IS a signal equal to the change in conduction angle proportional to the desired change in 
TCSC reactance to maintain zero error input The TCSC controller equation is gi’ven bv 

Xi — K[ Ky — Ki Iref 

The TCSC voltage can be expressed by its D — Q components as 


Ypcsc 



Figure 4 3 TCSC Controller 


TCSC TCSCD ' TCSCQ 


Linearizing the above equation we get 

(u 


“ i'^TCSCDO + ^'^TCSCd)'^ + i'^TCSCQO + ^'^TCSCq) 


neglecting the higher order terms we can write 


Au 


TCSC 


^TCSCDO 


'^TCSCO 


An, 


TCSCD 


'^Tcscqo A 
- f * ^^tcscq 

'^TCSCC 
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Similarly we can obtain the incremental relationship for line current The linearized equa 
tion for the line current is given bj 


= ^ 22 = A^aD + ^ A^oQ 

I'GO ^GO 

Now we can write the linearized TCSC control system dynamic equations in the state space 
form as 


Axc — [Ag Axg + [Be] Aue 


(4 8) 


where, 


Ag ] = [0] and [ Be ] = 


Ki K K[ K v^QSQQ 

“^TCSCO ^TCSCO 


Kj Z<7D tQQ 

IGO 


Axg = [Aaii] and Aug = Aigd Aigq]^ The output equation is 

given by 


Aye = [Cg] Axg + [Dg] Aug 


(4 9) 


where Aye = Act The matrices [Cg] and [Dg] are given by [Cg] = [1] and 

Kp K Kp K Vtcscqo Kp inno '■gqo 

^TCSCO ’^TCSCO 

4 13 Combined Generator - Network System Model 

The overall system mcludes subsystems of generator, network and TCSC controller The 
interconnection pattern of these subsystems is shown in Fig 4 4 The generator system 
input quantities are obtained from the network output variables, and the network input 
quantities are obtamed from the generator system output variables and the TCSC con 
troller output variables The TCSC controller inputs are related to generator and network 
outputs These subs> stems are interconnected as shown m Fig 4 4 to obtain overall system 


[Dc\ = 
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model The state and output equations of the total system can be written in the following 
state space form 


Ax = [A] Ax + [B] Au 

(4 10) 

Ay = [C] Ax + [D] Au 

(4 11) 


where Ax = [Axg Ax^ Axg]^, An = [Au^ Au^ Au^]^ Ay = [Ayg Ay^ Ay^f 
and the matrices 


Ag 


■ Bg 

Av 

[B] = 

By 

Ag _ 

Be 


■ Cg 


"Dg 

Cjv 

and [ D ] = 


Cg . 

Dg 


For the study system the input and output vectors of each subsystem are gi^ en by 
Aug = [AvgD Avgof Aya — \_Aigd Aiqq Aiqd ^‘^gq ] , 

Auiv = [Aig£) Aiqq Act Aigd > ^Yn = [Au^xj Au^q Avccd ^'^^ccqY' , 

Aiix7 A^Vqqq CLTid jAy Q — A^g 

The input to each subsystem is related to the output variables of the other subsystems 
The equations for various subsystem inputs can be combined to give the overall system 
input vector The relationship between the overall system input vector and the output 
vector IS expressed through interface matrix as 

Au = [F] Ay (4 12) 

Now utilizing eq (4 12) and then combining eqs (4 10) and (4 11) with the mterface 
matrix [F] we get the overall system model as 

Ax = [A + BF(I - DF)"^C] Ax = [At] Ax (4 13) 
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The details of the [F] matnx are gi'ven m Appendix E 


Generator system AC Network 



Figure 4 4 Interconnection Diagram 


4 2 Eigenvalue Analysis 

The study of lEEE-FBM is earned out to examine the performance of the controller m 
damping SSR oscillations through real power modulation The study system configuration 
IS shown in Fig 4 1 and the data is given in Appendnx B The generator is delnenng 0 9 
per unit power to transmission system at 0 9 power factor lag The field voltage is assumed 
to be constant and governor action is not considered Self damping of 0 2 per umt and 
mutual damping of 0 3 per unit is added to the shaft [13] The target compensation of 
the line is assumed to be 70% (0 35 per unit) This is achieved through a fixed series 
compensation (Cp as shown m Fig 4 2) which is assumed to be of 0 2546 per unit ( 72 74% 
of target compensation) and the b dance can be (0 0954 per unit) achieved through TCSC 
which has vernier control (enabled through a combination of fixed capacitor and a mductor 
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connected m senes with anti parallel thyristors) The ratio of fixed capacitor reactance 
and inductor reactance of TCSC module is chosen as 4 [39] 

The overall system eigenvalues are listed in Table 41 In this table trv o different sets of 
eigenvalues are gnen One with the entire compensation of 70 % being provided ^uth the 
fixed capacitor and the other parth through fixed compensation and the balance through 
the TCSC as shown m Fig 4 2 The proportional gam Kp = 0 01 and mtegral gam 
Kf = 10 are selected for the stud\ It can be observed from the Table 4 1 that mode 2 
is negatively damped for with entire fixed capacitance With the introduction of TCSC 
mode 2 is positively damped Thus TCSC contributes towards the damping of mode 2 
There is slight improvement m the mode 3 damping The damping of mode 1 is increased, 
whereas the mode 0 damping is decreased The subsynchronous network mode frequency 
IS increased, whereas the supersynchronous network mode frequency is decreased The 
damping of subsvnchronous network mode is increased 


Table 4 1 Eigenvalues of TCSC with continuous time control 


With Fixed 
Capacitance 

With TCSC 

Const Angle control 

(Ap = 0 01 Wj = 10) 

Comments 

4 41 ± j6l2 432 

4 71 ± j497 92 

Sups\Ti Net Mode 

3 44 ± jl41 11 

3 29 ± j255 93 

Subs\’n Net Mode 

1 38 ± j298 158 

1 38 ± j298 16 

Mode 5 

0 27 ± j203 03 

0 27 ± j203 24 

Mode 4 

0 8860 ± jl59 11 

0 92 ± jl59 48 

klode 3 

0 03 ± jl29 142 

0 042 ± jl28 88 

Mode 2 

0 10 ± j98 64 

0 151 ± j98 04 

Mode 1 

0 431 ± jl0 55 

3 91 

1 019 

0 21 ± j7 58 

0 ± j377 

3 18 +j 0 010 

0 04 + j 0 30 

0 08 j 0 30 

Mode 0 
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The effect of controller gams on torsional modes is shoum in Fig 4 5 From the Fig 
4 5, it can be observed that the controller gam Kj (with Kp fixed at 0 01) has minimum 
influence on torsional mode damping, whereas, the network mode damping is significantly 
influenced by the gam The damping of the network mode is decreased with the increase in 
controller gam Kj With Kj fixed at 10, the variation in controller gain Kp also has only 
minimal influence on the torsional modes However the network mode gets destabilized 
at a value of iifp = 0 095 



Figure 4 5 Variation of real part of system modes with Controller gams 

The effect of \ariation m TCSG reactance, which leads to variation in power flow 
over the line, is shown m Table 4 2 with Kj = 10 and Kp = 0 01 The damping of 
the subsynchronous network mode is increased with the increase in the value of TCSG 
reactance, whereas, the damping of supersynchronous network mode is decreased Also 
the damping of the zeroth mode is increased with the increase in the TCSG reactance 
From the eigenvalue analysis study, it can be observed that the TGSG with the constant 
angle control as suggested through the control of line voltage drop having more influence 
on the network mode damping This is due to the fact that the feedback signal considered 
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IS the electrical system quantity and it is expected that it influences the electrical modes 
compared to the mechanical system modes [13] 


Table 4 2 Eigenvalues for different TCSC reactances 


^Tcsc = 0 0932p u 
P=0 8 p u 

^Tcsc = 0 0954p u 

P = 0 9 p u 

^Tcsc — 0 1004p u 
P=1 0 p u 

Comments 

4 79 ± j497 94 

4 71 ± j497 92 

4 67± j497 89 

Supsyn Net Mode 

3 18 ± j255 8b 

3 29 ± j255 93 

3 5o± j255 99 

Subsyn Net Mode 

1 38 ± j298 168 

1 38 ± j298 16 

1 38± j298 16 

Mode 5 

0 27 ± j203 24 

0 27 ± j203 24 

0 27± j203 23 

Mode 4 

0 92 ±jl59 47 

1 0 92 ± jlo9 48 

0 92 ± j 150 47 

Mode 3 

0 042 ±jl28 88 

0 042 ±j 128 88 

0 042 ± j 128 88 

Mode 2 

0 15±j98 04 , 

0 151±j98 04 

0 lo ± j 98 04 

Mode 1 

0 201 ±j7o8 

0± j377 1 

1 0 21±j7 58 

‘ 0 ± j377 

0 23 ± J 7 58 

0 ± j377 

Mode 0 


4 3 Time Domain Simulation 

To validate the obser\ations made through the results shown m Fig 4 b that the controller 
gains have minimal influence on the torsional modes, a time domain simulation of the stud'v 
system has been carried out for two different sets of controller parameters A disturbance of 
0 1 per unit is applied at the HP turbine input for a period of 0 1 sec The variation in shaft 
torque and line power is shown in Figs 4 6 and 4 7 As expected the same real power flow 
on the line is achieved tvith two different sets of controller parameters The various shaft 
torques are almost identical with the two sets of controller parameters, thus, confirming the 
results shown m Fig 4 o 

The performance of the TCSC with the constant angle control which is achieved by 
maintaining the voltage drop across the line considered, is further examined with the 
syste n subjected to a major disturbance The disturbance considered is a three phase to 
ground fault of 4 5 cj cle duration at infinite bus end with a fault impedance of 0 04 per 


75 




unit The plots of the shaft torques are shown in Fig 4 8 It can be noticed from the 
Fig 4 8 that the response curves have decreasing amplitudes The same behavior can be 
observed from the generator real power response and the load angle response shown in Fig 
49 

The above results indicates that in the event of large disturbance the study system 
IS stable with the TCSC and that the system behavior is not adverseh affected by SSR 
From the response curves it can be observed that the peak values of Gen Exc torque 
LPB Gen torque and LP^-LPB torque are 0 315 per unit 2 25 per unit and 2 012 per unit 
respectively 

4 4 Conclusions 

In this chapter, the effect of TCSC on subs\Tichronous resonance phenomenon has been 
investigated The control of TCSC based on constant angle control strategy has been 
considered This control is achieved by mamtaimng the voltage drop across the line is 
constant It has been shown that the TCSC has a beneficial influence on the SSR response 
of the system while the TCSC controller parameters influence the network modes They 
have minimal influence on mechanical system modes The study results obtained through 
eigenvalue analysis are validated through detailed digital simulation 

For the purpose of eigenvalue analysis the TCSC compensated studv system has been 
modeled m the continuous time linearized domain and however, does not consider the 
discrete nature of TCSC operation The current trend seems to be explormg the possibility 
of application of digital control schemes for power transmission systems and this will be 
discussed m the next chapter 
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Figure 4 6 System response with Ki =20 and Kp=0 001 
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LPA LPB Torque in per unit 



Figure 4 7 System response with Kj =50 and Kp=Q 067 
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Figure 4 9 System response with Constant angle control based TCSC 
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Chapter 5 


SSR MITIGATION USING 
DIGITAL CONTROL OF TCSC 


In the previous chapter, the efficacy of comentional TCSC control has been examined in 
damping SSR oscillations There has howe\er been a trend to explore the suitability and 
use of digital control in controllable elements being employed in power transmission 

In this chapter, a novel digital control scheme is proposed for TCSC Generally, in TCSC 
control the approach is to utilize a look up table in place of closed form transcendental 
expression to determine the firing angle for the desired TCSC reactance This is preferred as 
the solution of the closed form expression in real time is both computationally tedious and 
time consuming The look up table is generated off line for a given TCSC In this chapter 
use of an estimated characteristic of the TCSC is suggested to relate the TCSC reactance 
to the firing angle The estimated characteristic is incorporated m a D — Q domain state 
space model of TCSC compensated transrmssion line The resulting nonlinear model is 
utilized to derive a linear discrete time domain state space model and the controller is 
designed using state feedback approach 

The discrete time domain model of the TCSC compensated transmission line along 
with controller is converted into a continue us time domain model to mterface with the 
synchronous machine model and obtain hemogeneous state space equation Using the 
system model so derived, the eigemalue anaJ^sis is carried out to observe the influence of 
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digital controller on torsional mode damping The results obtained through the eigenvalue 
analysis are validated with time domain simulation 

5 1 Estimation of TCSC Characteristic 

The TCSC consists of a capacitor in parallel with an mductor connected to a pair of anti 
parallel thyristors By varying the firing angle of the thvristors the inductor reactance 
IS varied and which in turn changes the effective reactance of the TCSC Thus, we can 
represent the TCSC as a variable impedance device The TCSC is generally designed to 
operate in the capacitive region, while inductive mode operation can be used during the 
fault conditions In this section, the estimation of the TCSC characteristic in the capacitive 
operating region is described 

Let us denote the TCSC capacitor reactance by Xc the thyristor controlled inductor 
reactance by Xp and the thyristor controlled reactor filing angle by a The overall TCSC 
reactance {Xpcsc) is given by [22] 

^TCSC ~ ~ l?l(-^C + 172) + ^4175 


where 


^1 

m 

Vi 

Vh 


2 ( 7 r - oi) + s%n2{r — a) 

TT 

XcXp 
Xc — Xp 



Tj^tan {773 (tt - q;)} — tan(7r - a) 

4r]2Cos'^{% - a) 
irXp 


T1 e variation in TCSC reactance with a is shown m Fig 5 1 This curve can now be 
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approximated by an exponential function which is given by 


Of 

^Tcsc = + ^0 exp( ) 

Ofc 


(d 2) 


where 


Xj'QsCmax 

= p ■ 

p = exp(-^) 

Xroscmax IS the maximum value of Xtcsc ^mm is the corresponding firing angle and 
ac IS called the a constant of this device The a-constant is defined by the TCR firing 
angle a which is required to achie\e the TCSC reactance of Xa + 0 632{XTCscmax — ^c) 
Consequently, etc is obtained as 0 632 The estimated curve is also shown in Fig 5 1 


5 2 TCSC Compensated Transmission Line Model 

The system under study is shown in Fig 5 2 This is essentially the IEEE FBM with 
a fixed series capacitor and a TCSC connected to it Let us denote the reactance of the 
fixed capacitor and TCSC respectively by Xp and Xtcsc The system equations are then 
written in the synchronously rotating D — Q axis reference frame as 

[x] = [F] X + [G] V 3) 

where the state vector x contains the D - Q axes line currents and capacitor voltages i e 
x^ = [zii, iLQ vcD 'iJeg] and the vector v is given by v^ = [vgo - v^d '^gQ ' ^ooq] and 
the matrices [F] and [G] are given by 

0 i 
0 0 
0 0 _ 

It ,s to be noted that X,. = 3 ^, Xrese = 25^;? and L is the total inductance of the 
transmission line and the transformers Combining eqs (4 2) and (4 3) we get 

X = Ax + Bv + Cxu 


[F] 


-5 

L 

UJ 

Cp ‘ Ctcsc 

n 


“O; 

R 
" L 


+ 


/■y 


_-i 

L 

0 

0 

u 


0 

1 

L 

—U) 
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where [B] [G], u exp( ^ ) and the matrices [A], [B] and [C] are given by 


[A] = 


R 

L 


U 


-W 
R 
' L 


U){Xp + Xcmin) 0 

0 ‘^(^F+Xcmm) 


1 

L 

0 ■ 


■ 1 

L 

0 ■ 

0 

0 

1 

L 

—a; 

.1B1 = 

0 

0 

1 

L 

0 

u 

0 


0 

0 _ 


[C] = 


0 0 0 0 

0 0 0 0 

^beta^ 0 0 0 

0 Ki^eta^ 0 0 


Since u IS d. scalar, the system eq (5 4) reduces to 

X = ^A + Cexp(--^'j x+Bv 


(5 5) 


The matrix 


[A + C exp{—-^) has two pairs of eigen\-alues It is observed from the root 
locus plot of the eigenvalues of the above matrix for vanous values of a, that the real parts 
of the two pairs of eigenvalues are constant as a changes from 180° to 140° and one pair 
of eigenvalues moves towards the real axis while the other pair moves outwards [40] 


5 2 1 Linearized Model of TCSC Compensated Transmission 
Line 


The solution of eq (5 5) is given by 

x{tj) = exp{(A + Cu)AT}x(to) exp{(A+ Cu)(t/ - t)}B v(r)(ir 
= exp (FAT) x(to) + exp {r(t/ — r)} B v{T)dT 


(5 6) 


where F = A + C exp j > = tf — to, and to ond tf are initial and final values of 

the sampling instants Equation (4 6) is nonlinear due to the presence of the term u that 
IS a function of the firing angle We now have to linearize eq (5 6) To do that we note 
that th independent variables are x(to), v and a Expanding eq (5 6) in Taylors series 
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around nominal values of xq vq, and neglecting the second and higher order terms -we 
get 


(5 7) 

= AAx(to) + B Av(to) + CAa{to) 
where Ax Av and Aa denote the perturbations from the respective nominal values The 
derivation of the matrices A B and C are given below From the eq (5 6) we notice that 
the second term is independent of x{to) Thus 



= A == AA = exp(rAr) 

uyi{to) 

Similarly the first term on the right hand side of eq (5 6) is independent of v Thus, 

= B ~ [ exp {r(AT — r)} Bdr 

Jo 

Since r IS non singular, the above equation can be simplified as 


B = AB =r-n^A-J)B 


Now the matrix C is given by 


da 


C can be divided into two parts Ci and C 2 such that C — AC — Ci + C 2 These two 
terms are given below 


Cl 


A 

da 


a 


exp 


A + C exp( ) / x(to) 


ar 


AT. , ao. 

AACexp( ) Xo 

Oc Gic 


r\ pAT 

C^ = A exp{A + Cu)(t/-r)}B v(r)dr 
uOi J 0 

Assuming v(t) = Vq, 0 < ^ < AT and solving we get 

C 2 = {ATexp(rAT) - exp(rAT) + C B 

etc 
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It IS assumed that the firing angle is generated based on the individual firing angle control 
such that a is updated six times in each cycle Now representing the eq (5 7) in discrete 
time domain as 


£\x{k + 1) = AA Ax{k) + AB Av(;c) 4- AC Aa(k) (5 8) 

It is to be noted that since the infinite bus voltage is assumed to be constant we can write 

Av^ = [AUgl, AUgg] 

5 3 TCSC Control 

The discrete time model of the system is used with a pole placement regulatmg control 
law of the form 


Aaik) = - [K] Ax{k) (6 9) 

where [K] is a feed back gam matrix chosen using the pole shift control technique [41] In 
this method the open loop system poles are radiallv shifted to more stable locations in 
the closed loop Let the open loop characteristic equation of a discrete time s\'stem be 
given by 

Gol{z) = + 52-2"“^ H 9n-iz + 

then in the pole shift control, the closed loop characteristic equation is specified by 

Gciiz) = z" + 51 A 0"-^ + 52 A^ 2:"-^ + + 5n-i A”-^ + 5n A" 

where A (0 < A < 1) is the pole shift factor It is to be noted that as A reduces, the penalty 
on the control also reduces 

The TCSC control system block diagram is { iven in Fig 5 3 The three phase line 
currents and capacitor voltages are transformed 1 1 the D — Q domain to obtain the vector 
X The vector x is compared with the nominal value to obtain Ax The output Ao; 
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of the controller is added to the nominal value to obtain TCSC firing angle \ow the 
discrete time TCSC model incorporating the controller can be written as 

Ax(A: 4- 1) = (AA - ACK) Ax(yt) + AB Av{k) (5 10) 

To interface the above discrete time model with the continuous time turbine generator 



Figure 5 2 The schematic representation of TCSC control system 


system model for eigenvalue analysis, we convert the above equation to continuous time 
and obtain the following equation 


Ax = 



Ax + 



Av 


(5 11) 


where 


■ A' 

B' ■ 

= — ^ log 

O 24 

O 22 

AT ^ 


AA-ACK AB 
O2 4 I2 2 
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In the above equation 0, j is a null matrix with % rows and j columns and I 2 2 is an identity 
matrix of order 2 It is to be noted that the sampling time AT defined before is used for 
the above transform 


5 4 Combined Generator-Network System Model 


The basic purpose of the generator network interface is to preserve the structure of the 
network given by eq (5 11), since this will allow us to check the validity of the state 
feedback controller vis a vis the stability of the system To facilitate such an interface we 
define a fictitious capacitor of value Cn at the generator terminals The value of Cn can be 
chosen to be arbitrarily small (Ac = 100 per unit) as this is only used for interfacing The 
network equations then can be written from eq (5 11) as 


Axiv = [Aw] Axw + [Bw] Auw 
= [Cw] Axw 


(5 12) 


where Axw = Iyiiq I^vcd ^'^cq Augi? Augg]^, Auw — [Aigx) Aigg] and Ayw 


[Augip iYVgqf and the matrices are 


[A/^] 


A 

0 


[Ai] 


0 - 




B 

Ai 

0 0 0 0 
0 0 0 0 


[Bw] = 


0 

0 

0 

0 

Cn 

0 


0 

0 

0 

0 

0 

i- 

Cn 


[Cw] 


0 0 0 0 1 0 
0 0 0 0 0 1 
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The state and output equations of generator system are mentioned in Chapter 2 (Section 
2 1) and are reproduced here 


Axg = [Ag] Axg + [Bg] Aug 

(5 13) 

Aye = [Cg] Axg 

where Axq — [Ax^j^ Aa;^], Aug = [Auj£) Avgof , Ajq = [Aid Aiq]^ Axjx is the state 
vector of stator and rotor variables and Ax^ is the state vector of mechanical svstem 
variables It is assumed that machine output currents Aid and Aiq are equal to the 
network injection currents i e , Aiqd and Azqq respectively It is to be mentioned here 
that the above output equation of the generator system is different to the one given in 
Chapter 2 (Section 2 1) Accordingly the dimension of the [Cg] matrix will change and 
also [Cg] = [Cgi] The matrices [Ag], [Bg] and [Cgi] are defined in Chapter 2 We now 
combine eqs (4 12) and (4 13) to obtain the overall svstem model as 

Ax = [At] Ax (5 14) 

where Ax^ = [Axg Axg] and 


. ] 


1 - 

Ag 

1 

“Bg 



J - 

Ajv _ 






5 5 Eigenvalue Analysis 


The studv of a lEEE-PBM is carried out to examine the performance of the proposed 
digital control in damping the SSR oscillations through real power modulation The study 
system configuration is shown in Fig 5 2 and the data is given in Appendix B The 
generator field voltage is assumed to be constant Series capacitive compensation is 70 9% 
of transmission line reactance Out of this 40% is in the form of fixed capacitive reactance 
Cp and the rest is in the form of TCSC The generator is deh vermg 1 43 per unit power to 
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the infinite bus Self damping of 0 2 per unit and mutual damping of 0 3 per unit is added 
to the shaft to account for steam damping [13] 

A pole shift controller is used for the state feedback control of eq (5 9) In this method 
the open loop system poles are radially shifted to more stable locations in the closed loop 
The pole shift factor value can be selected m the range between 0 and 1 As mentioned 
previously in Section 5 3, the control effort is less when the value of pole shift factor is close 
to 1 and the reduction in pole shift factor ■value increases the control effort The chosen 
value of pole shift factor must result in the stable system operation 

System eigenvalues are listed m Table o 1 for pole shift factors of 0 95 and 0 8 It 
IS observed that the damping of supersynchronous and subsynchronous network modes 
are increased with the introduction of TCSC as compared to that obtained with fixed 
compensation For pole shift factor of 0 95, mode 1 and mode 2 are unstable All torsional 
modes are damped with a pole shift factor of 0 8 With the decrease of pole shift factor to 
0 8, the mode 1 and mode 2 are damped but the damping of mode 0 is slightly reduced 
Also, the damping of supersynchronous and subsvnchronous network modes is considerably 
increased The frequency of oscillation of subsvnchronous network mode is increased with 
decrease of pole shift factor value from 0 9o to 0 8 The subsynchronous network mode for 
a pole shift factor of 0 8 has a time constant of 11 6 ms whereas without TCSC the value 
IS around 0 5 sec This indicates fast damping of subsvnchronous network mode 

The requirements for SSR mitigation are either to mcrease mechanical damping, or 
reduce the electrical destabilization One of the major concerns of SSR control is the res- 
onant response of the torsional system m the post fault state In the case of undischarged 
fixed capacitor, while returning to the post fault equihbnum state, it produces subsyn- 
chronous network currents and associated generated electromagnetic torque that may have 
frequ ncy components close to torsional frequency Even if the post fault state is stable, 
the elfect of resonance causes very large shaft torques to build up The rapid discharge 
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Table 5 1 Eigenvalues of TCSC with digital control 


Fixed 

capacitance 

Pole shift 
factor of 0 95 

Pole shift 
factor of 0 8 

Comments 

4 41 ± j612 432 

22 79 ± j613 26 

85 69 ±j612 13 

Supersyn Net Mode 

3 44 ±jl4111 

22 10 ± jl37 74 

81 69 ± 3145 74 

Subsyn Net Mode 

1 380 ± j298 15 

1 380 ± j298 15 

1 380 ± j298 15 

Mode 5 

0 27 ± j203 03 

0 244 ± j203 027 

0 242 ± j203 10 

Mode 4 

0 886 ± j1d9 11 

0 787 ± j159 29 

0 88 ± jl59 43 

Mode 3 

0 03 ±jl29 142 

0 079 ± jl28 92 

0 023 ± jl28 875 

Mode 2 

0 10 ± j98 64 

0 088 ± j98 46 

0 042 ± j98 04 

Mode 1 

0 431 ± jlO 00 

0 3024 ± j9 938 

0 2618±j8 16 

Mode 0 


2 93±jl798 95 

4 67± 3 1797 64 



3 57 ± jl378 81 

2 93 ± 31378 47 


3 91 

40 

3 84 


0 19 

0 11 

0 058 



of capacitor energy by TCSC can eliminate subsynchronous currents quickly The time 
constant of subsynchronous network mode for a pole shift factor of 0 8 is less than one 
cycle This indicates that the subsynchronous currents are damped quickly Also the in- 
teraction between the torques produced by the subsynchronous currents and shaft torques 
IS minimum 

5 6 Time Domain Simulation 

To obser-e the effectiveness of proposed TCSC controller, time domain simulation of the 
study system is earned out using PSCAD/EMTDC package version 3 0, by includmg the 
system non linearities The field voltage is assumed to be constant A 4 5 cycle, three 
phase to ground fault with 0 04 per unit fault impedance at the infinite bus is applied 
when the system is in steady state From the eigenvalue analysis, it can be observed that 
mode 1 and mode 2 are negatively damped for pole shift factor of 0 95 Fig 5 3 shows the 
shaft torques when TCSC is operating in feedback control mode with a pole shift factor of 
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0 95 It can be seen that the shaft torque oscillations are increasing These results are in 
agreement with the eigenvalue analysis results 

The plots of the shaft torques with pole shift factor of 0 8 are shown m Fig 5 4 It can 
be observed from Fig 5 4 that with the introduction of TCSC all shaft modes are damped 
The settling time of LPA LPB torque is relatively large because of low damping The 
same observation is evident from the eigenvalues shown in Table o 1 where the damping 
of mode 2 (which corresponds to LPA LPB torque) is poor The plots of generator real 
power and load angle are shown in Fig 5 5 for this stable case 

5 7 Conclusions 

A novel digital control scheme is proposed for TCSC control To study the performance of 
the control scheme a discrete time model of the TCSC has been developed and which is 
suitablv interfaced with the generator system model The control is based on state feedback 
approach which utilizes capacitor voltage and current Eigenvalue anah sis is carried out 
to examine the performance of the proposed digital control with TCSC in damping the 
SSR oscillations The variation in pole shift factor affects the damping of mode 1, mode 
2 supersynchronous and subsynchronous network modes The TCSC contributes to the 
damping of network modes to a great extent, but at the same time it also reduces the 
damping of some of the stable torsional modes marginally The choice of pole shift factor, 
therefore, becomes critical Time domain simulations are carried out to study the transient 
torque case It is shown that the introduction of TCSC results in the stabilization of 
torsional modes with proper choice of pole shift factor 
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Figure 5 3 System response with digital control based TCSC ( pole shift factor of 0 95) 
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Chapter 6 


TORSIONAL INTERACTION 
BETWEEN TCSC AND PSS 


In the study presented in Chapters 4 and 5, the effect of excitation system and power 
system stabilizer is not considered The power system stabilizer (PSS) is generally used to 
damp the low frequency oscillations with generator speed as auxiliary signal The generator 
speed has the components of shaft modes and these components have significant effect on 
the performance of power system stabilizer When PSS is used with generator speed as 
supplementary signal there may be a possibility of it interacting with TCSC control m the 
range of torsional mode frequencies The control interaction studj is, therefore, important 
to ensure proper design and coordination between various controllers 

Power system stabilizers (PSS) can extend stability hmits of systems that exhibit lightly 
damped oscillations in the low frequency (0 2 to 2 5 Hz) range This is accomplished 
via excitation control Torsional oscillations of turbine generators are inherently lightly 
damped There have been several instances of torsional mode instability due to interaction 
between generatmg unit excitation and nearby HVDC converter controls [42] In the 
previous chapters it has been shown that the TCSC with suitable control strategy provides 
positive damping to unstable torsional modes However, presence of PSS may influence 
this performance In this chapter, the extent and nature of control mteraction between PSS 
and TCSC has been studied The TCSC considered to be equipped with the controllers as 
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described m Chapters 4 and 5 


6 1 System Representation 

As before, the study system considered is IEEE FBM For this system the various com 
ponent models viz generator system model network model and TCSG model have been 
developed m the previous chapters These models are utilized here to develop the complete 
system model while suitably incorporating the excitation system and PSS 

6 11 Excitation System and PSS Model 

Consider a static exciter that includes a single time constant voltage regulator as shown 
in Fig 6 1 [4] In the block diagram of the exciter Vg is the machine terminal voltage 
and Vs is the output signal from auxiliary controllers such as PSS The excitation system 
dynamics can be written as 

EFD = ~E^D + ^{v„f-v, + v,) (61) 

-‘■R J-R 

where Tn is the exciter time constant Also Vg is related to its D and Q axis components 
by 

^9 = (b 2) 

Substituting eq (6 2) in eq (6 1) and linearizing we get 

ASfo = + ^A», - + ^Av„/ (6 3) 

The PSS consists of a washout circuit, lead/lag blocks and a limiter Its block diagram 
IS shown in Fig 6 2 The input signal to the PSS is the generator shp For the sake 
of convenience in writing the state equations for PSS the block diagram of Fig 6 2 can 
be modified to the form shown in Fig 6 3 From this figure, it can be noticed that the 
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Figure 6 1 Static Exciter 

system is having three poles and requires three state variables to descnbe its dynamics 

The following linearized state equations for the PSS can therefore be written as 

Axi = —GiAxi + C4AX2 + KsCs^x^ + KsCs^Sg 

i\X2 — — C\/\X2 d" d" KsC^iC^Sg 

11 

Axs = - T^/^Sg 

Iw 

~ ^1 d" C 2 X 2 d" K 3 C 2 XZ d" Ks^2^3 

where Ci = ^ C2 = ^ C3 = 1 - C4 = C1C3 C5 = CiCzCz and T;v is the washout 
time constant Combining eqs 6 3 and 6 4 , we can write the following state space equation 
for the static exciter with PSS 

Axa = [Aei] Axei d- [Ben] Auei d- [Bei 2 ] Aug + [Eei] Av„/ (6 5) 

Ayel = [Cel] Axei (6 6) 

where Ax^ = [A^^^n Axx Axz Ais]^, Auei = As^ and Ayei = /^Efd Smce Vref 
IS assumed constant, ^Vref — 0 matrices [Aei], [Beii]) [Bei 2 ], [Eei] and [Cgi] are 

defined in Appendix F 
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0 

ST„ 


L(1 + STi)' 

! 

Vs 


1 + ST, 

“ 

! 

(1 + stD' i 


' > 


Washout Lead/Lag Limiter 

Figure 6 2 Block diagrELm of Power System Stabilizer 


6 12 Complete Generator System Model with Excitation system 

Defining state -vector as x^J, and combining eqs (2 19) (6 5) and (6 6), 

we can -write the following state space equation 



Asi 


Ball 

Axg = 


Axg + 

0 


Agi 


0 


Ausu+ 


Bal2 


■ Bal3 ■ 


0 


0 

0 

0 

Au^12 -+- 

0 

. . 

Aug + 

I 

Au„i + 

0 

. Bell _ 


AUei 


(6 7) 


Note that = ATe = Cmei Ax^i, AUsii = Aym = C^AXm, AU512 — Ayei — CeiAXei, 
AUel = Cr„2Ax^ and Cr „2 = [0 1] Hence we can rewrite the above equation as 


Ax<3 = 


Asi 

BsIiCtu 

Bai2Cei 


Bai3 

BM^mel 

Am 

0 

Axg + 

0 

0 

BeiiCm2 

A-el 


. Bel2 . 


A-Ug 


= [ Ag ] Axg + [ Bg ] Aug 

where Aug = [Aigo Avgqf The output equation of the generator as defined m Chapter 
2, IS 

AyG = [Cg] Axg + [Dg] Aug (6 9 ) 

It may be noted that the dimensions of matrices [Cg] and [Dg] w 11 get modified suitably 
in view of the inclusion of the excitation system and PSS model 
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6 2 Combined Generator - Network System Model 

Constant angle control based TCSC 

State space model of TCSC with constant angle control is derived in Chapter 4 (Section 
4 1) This IS reproduced below 

Axc = [Ac] Axc + [Ba] Auc 
Aye = [Cc] Axc 

where Axc = [Axi AX 2 ] , Auc = [^vtcscd ^'^tcscq ^t-gd Azgq]^ and Aye = Act 
T he network model as derived in Chapter 4 is 

Axiv = [Aiv] Ax;^ 4 - [Bat] Aujy 
Ayw = [Civ] Axiv + [Cat] Aua^ 

where 


Axat — [Avcfd ^'^cfq ^vtcscd ^VTcscof , Auat = [Alee Aicq Act Aiee Aieg]^ 

and AyA/- = [Augo Au^g ^vtcscd ^vtcscqT Following the usual procedure network 
and TCSC system equations can be combined with eqs (6 8) and (6 9) to get the overall 
system model in the form Ax = [At] Ax 
D igital Control based TCSC 

The network incorporating the model of TCSC based on digital control strategy has been 
derived in Chapter 5 (Section 5 4) For the sake of completeness the network model is 
given below 

Axa/- = [Aw] Axat + [Bw] Auat 
AyAT = [Cm] Axw 


where 


Axw = [Ai£,e Aix,g ^vcd ^Vgn Aupg]^ 
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Aujv - [A 2 G 0 A^GQ]^ and = [Avgo Au^q]^ 


Following the usual procedure the network equations can be combined with eqs (6 8) and 
(6 9) to get the overall system model in the form Ax = [At] Ax 

It may be noted that the overall system matrix [Ay] used abo\e has the same structure 
as the matrix [Ax] given in Chapters 4 and 5 except that its component matrix [A^] now 
gets modified as per eqn (6 8) to consider inclusion of excitation svstem and PSS 

6 3 Case Studies 

To investigate the effect of excitation system (Automatic Voltage Regulator (AVR) and 
PSS), the following three situations are considered with the study system augmented to 
include 

(1) Fixed series compensation 

(2) Constant angle control based TCSC 

(3) Digital control based TCSC 

The study system is essentially the IEEE FBM The study is cairned out using eigenvalue 
analysis The excitation system and PSS data is given m Appendix F and the FBM data 
is given in Appendix B 

6 3 1 Fixed Series Compensation 

Before proceeding with the analysis of the effect of AVR and PSS with TCSC controls, 
it would be beneficial to first examine the generator excitation controls with fixed series 
compensation Considering fixed series compensation in the study system the following 
three scenarios are examined 

(i) without AVR and PSS 

(ii) with AVR and without PSS 
(ill) with AVR and PSS 
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The eigenvalue analysis result for these cases is given in Table 6 1 In the present study 
the stabilizer gam is chosen as K 3 = Z 

There are five pairs of torsional modes which are important from the view pomt of sta 
bilizer application one pair of each supersynchronous network mode and subsynchronous 
network mode and one pair of zeroth mode The sixth mode identified as exciter mode 
primarilv associated with field voltage and exist only in the presence of speed stabilization 
This mode is usually well damped like network modes 

From the eigenvalue analysis results given in Table 6 1 it can be seen that mode 2 
IS unstable without AVR and PSS (case 1 ) This mode becomes stable when the A.VR 
is introduced (case 11 ) However in this case the zeroth mode is destabilized The speed 
signal IS the obvious choice for PSS to damp the zeroth mode Further introduction of 
PSS (case 111 ) stabilizes this mode, while modes 1 and 2 become unstable Thus, the 
PSS has adverse effect on torsional mode damping while contributing to the damping of 
zeroth mode Furthermore, the damping of modes 3 and 4 is reduced The speed signal is 
inherently sensitive to the presence of torsional oscillations and this is leading to negative 
damping of some of the torsional modes 

From the results it is clear that the PSS has significant effect on torsional mode stabihty 
It IS to be mentioned here that at lower values of PSS gam (1 e Kg < 1) the torsional modes 
are not destabilized, but the damping of zeroth mode is poor The primary objects e of 
the PSS IS, however, to damp the zeroth mode Destabilization of torsional modes due to 
PSS IS also observed in [43], but in case of an uncompensated transmission system 

6 3 2 Constant angle control based TCSC 

The effect of the stabilizer on torsional modes m the presence of TCSC control is analyzed 
The integral gam Ki and proportional gam Kp of TCS 3 controller are set at values of 10 
and 0 01 respectively, and PSS gam is varied The res ilts are given m Table 6 2 It can 
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Table 6 1 System eigenvalues with fixed senes compensation 


Without AYR 
and PSS 

With AYR and 
without PSS 

With AYR 
and PSS 

comments 

4 42± j612 43 

3 45± jl41 11 

1 38 ± j298 16 

0 28 ± j203 03 

0 89 ± jl59 12 

0 03 ± jl29 14 

0 11 ± j98 651 

0 43 ± jlO 56 

3 91 

0 20 

4 4 ± j612 47 

2 72 ± jl41 64 

2 04 ± j298 16 

0 27 ± j203 03 

0 99 ± jl59 11 

0 019 ± jl29 13 

0 28 ± j98 60 

0 05 ± jlO 96 

23 35±jl9 03 

2 847 

4 406 ± j612 47 

2 76 ± jl41 63 

1 98 ± j298 16 

0 2 ± j203 03 

0 062 ± jlo9 12 

0 02 ± jl29 12 

0 08 ± j98 60 

0 043 ±jl0 62 
23 15±jl8 95 

32 99 

20 6 

0 87 

0 1 

Supersyn Net Mode 
Subsyn Net Mode 
Mode 5 

Mode 4 

Mode 3 

Mode 2 

Mode 1 

Mode 0 

Exciter mode 


Table 6 2 System Eigenvalues with Constant angle Control based TCSC 


PSS gam 

Ks=2 

PSS gam 

K,=-L 

PSS gam 

a:,=io 

comments 

4 857 ± j497 878 

3 301 ± j2oo 95 

1 38 ± j298 158 

0 24 ± j203 26 

0 90 ±jlo9 49 

0 03 ± jl28 892 

0 0795 ± j98 149 

0 316±j9 467 

21 62± j21 44 

0 ± j377 

2 8 +j 0 000 

0 005 + j 0 00 

0 10 jOO 

4 8563 ± j497 878 

3 379 ± j255 95 
138 ±j298 158 

0 22 ± j203 274 

0 88 ± jlo9 503 

0 0254 ± jl28 90 

0 005 ± j98 20 

1 133 ± j9 79 

20 71± j20 75 

0 ± j377 

2 8 +j 0 00 

0 005 + j 0 0 

01 J 00 

4 853 ± j497 88 

3 409 ± j255 91 

1 38 ± j298 158 

0 164 ±j203 304 

0 83 ±jl59 525 

0 003 ± J128 91 

0 26 ± j98 38 

4 66 ± jlO 91 

17 33±jl6 85 

0 ± j377 

2 75 +j 0 0 

0 005 + J 0 0 

0 1 J 0 0 

Supsvn Net Mode 
Subsyn Net Mode 
Mode 5 

Mode 4 

Ylode 3 

Ylode 2 

Ylode 1 

Mode 0 
Exciter mode 


be seen hat with an increase in PSS gam, the dampmg of modes 1 and 2 decreases For 
a PSS g in = 4, mode 1 has positive real part, indicating small signal instability while 
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both modes 1 and 2 become unstable for PSS gam of 10 Also, the damping of mode 3 and 
mode 4 IS decreased with the increase of PSS gam These modes however remains stable 
The damping of the zeroth mode increases with stabilizer gam At a PSS gam value of 
all torsional modes are stable It can be noted that with fixed series compensation the 
PSS has greater detrimental effect on the damping of modes 3 and 4 as compared to the 
present case In the presence of TCSC, PSS has detrimental effect on modes 1 and 2, but 
this effect is more pronounced at higher PSS gams 

6.3.3 Digital Control based TCSC 

The primary objective of the TCSC is to provide dampmg to torsional modes From the 
eigenvalue analysis results presented m Chapter 5, it is observed that the study system is 
stable for a pole shift factor of 0 8 Here the effect of PSS on torsional modes is exammed by 
varying PSS gams while keeping the pole shift factor at a value of 0 8 The results are given 
m Table 6 3 It can be seen that with an increase m PSS gam, the dampmg of modes 1 and 
2 decreases For a PSS gam K, = 2. mode 1 has positive real part, indicating small signal 
instability while both modes 1 and 2 become unstable for PSS gam of 4 Furthermore, the 
dampmg of mode 3 and mode 4 is decreased with the mcrease of PSS gam These modes 
however remains stable On the positive side, the dampmg of the zeroth mode mcreases 
with stabilizer gam At a PSS gain value of 0 5 all torsional modes are stable, but the 
dampmg of mode 1 and mode 2 is poor 

6.3.4 Discussion 

Prom the case studies presented m the previous sections the following general observations 
can be made 

1. The PSS has greater positive influence on zeroth mode(power swmg mode) m the 
presence of TCSC as compared to the fixed senes .compensated system. This is 
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Table 6.3: System Eigenvalues with Digital Control based TCSC 


PSS gam 

K,=4 


comments 


Supersyn Net Mode 


PSS gam 
is:, = 05 


PSS gam 

K, = 2 


-85 65± j612 37 
-80 624±jl45 32 
-1 38 ± j298 15 
-0 241 ± j203 12 
-0 8798 ± jl59 44 

-0 023 ±j 128 80 

-0 039 ± j98 08 
-0 0044 ± j9 06 
-23 41 ± jl9 55 
-3 79 ±jl798 19 
-2 22 ±jl377 30 
-33.64 
-29 064 
-2 844 
-0 1 


-85 65± j612 37 
-80 68± jl45 36 
-1 38 ± j298 15 
-0 225 ± j203 11 
-0 861 ± jl59 44 
-0 on ±jl28 88 
0 063 ± j98 13 
-0 8 ± j9 466 
-22 47±jl8 82 
-3 79 ±jl798 19 
-2 2 ±jl377 3 
-36 36 
-26 88 
-2 78 
-01 


-85 64± j612 38 
-80 75± jl45 41 
-1 38 ± j298 15 
-0 20± j203 11 
-0 83 ± jl59 43 
0 005 ± jl28 88 
0 1994 ± j98 21 
-2 0±jl0 06 
-21 16±jl7 73 
-3 79 ±jll798 19 
-2 2 ±jl377 3 
-38 88 
-24 83 
-2 7 
-0 1 


Subsyn Net Mode 
Mode 5 
Mode 4 
Mode 3 
Mode 2 
Mode 1 
Mode 0 
Exciter mode 


expected, as the basic function of the PSS is to damp the zeroth mode 

2 In the presence of TCSC, PSS has detrimental effect primarily on some torsional 
modes ( mode 1 and mode 2), but this effect is prominent at higher values of PSS 

gains 

To arcumveat the detrimental impact of PSS on torsional modes, it is possible to use [4] 
torsional filter m the speed signal loop of PSS This aspect is exammed m the next section. 

6.4 Excitation System with PSS and Torsional filter 

In this section the mathematical model of synchronous machine with excitation system, 
PSS and torsional filter is developed The schematic diagram of PSS mth torsional filter 

IS shown m Fig 6.4 The output of the PSS is denoted by V .i Prom this figure, it can 

be noticed that the torsional filter has two poles reqninng two state variables to describe 
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its dynamics and three state variables are associated with PSS as mentioned before The 
state space equations for the static exciter have been given earher m Section 6 1 The state 



Washout PSS Torsional filter Limiter 

Figure 6.4: Power System Stabihzer 

space equations for PSS with torsional filter are given by 

Xi = X2 

X2 = - 2CaJ„a;2 + (u^Usi 

0:3 = — C1X3 + CiX^ + KgCsXs + KsCsSg 

x^ = -CiX4 + K,C4X5 + KsC4Sg ( 610 ) 

I 1 1 

X5 = -7^X5-7f^Sg 

Iw -^w 

Vs = Xiul 

Vsi = 2:3 + C 2 Xi + KsC^xs + KgC^Sg 


where Ci = C2 = ft, C3 = 1 - ft, C4 = 0x^3, C5 = In the above equations 

the state variables aq and X2 are associated with torsional filter and xs, x^ and x^ are 
associated with PSS and washout circuit 

Now we linearize the above set of equations and combine with eq. (6 3) to get the 
following state space model of the excitation system that includes PSS and 1 orsional filter 

Axa = [Ael] AXel + [Bell] AUel + Peiz] Au^ + [Eel] AV„/ (6 H) 


108 




^y<i = [c,i] Ax,i 


(6 12 ) 

Where = [AE^o Ax, Ax^ Axa AX4 Ax^]^, Au,, = A5, and Ay., = AE,n The 
matrices [A.,], [B.,,], [Be,2], [C.,] and [E.,] are defined in Appendix F Since v^./ is 
assumed constant Avr./ = 0 Following the procedure given Section 6 1 we can write the 
state space equation for the combined generator system as 

Axg = [Ag] Axg + [Bg] Aug T 

AyG = [Cg] Axg + [Dg] Aug 

where Axq = [Ax^, Axm Ax.,]^, and Aug = [Axjg Au^q]^ Equation (6 13 ) can now 
be used to derive the complete system model as described m Section 6.1 considering both 
the constant angle control based TCSC and the digital control based TCSC 

6.5 System Study 

The overall system model as derived above is used to mvestigate the eflfect of torsional filter 
in the performance of PSS Two configurations of the study system axe considered One 
with constant angle control based TCSC and the other with digital control based TCSC 
The excitation system and PSS data is given in Appendix F and the FBM data is given 
in Appendix B 

6.5.1 System with Constant angle control based TCSC 

From the analysis presented in the previous section, it is observed that the PSS is desta- 
bilizing some of the stable torsional modes m the presence of TCSC control The system 
eigenvalues with variation m torsional filter parameters {ojn cind are given m Tables 6 4 
and 6 . 5 . Only torsional modes axe given in the table. It can be noticed from the results 
that all the torsional r lodes are stable In this study the PSS gam is selected as 10 How- 
ever, for this gam vah e mode 1 and mode 2 were unstable without torsional filter. This 
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shows that as expected, the filter makes the torsional modes insensitive to the speed signal 
stabilization 

The study is carried out to further examine the effect of filter parameters on the system 
performance The analysis considered the effect of variation of filter parameters (a;„ and Q 
Initiallv IS fixed at 0 5 and is varied For this three values of a!„ are considered varymg 
The results are given in Table 6 4 It can be observed from the results that the damping 
order of torsional modes is same for different values of Variation of Un has a sig nifi cant 
effect on zeroth mode and torsional filter mode damping However, its influence on the 
exciter mode damping is minimum Table 6 5 gives the results of variation of C with a;„ 
fixed at 7 rad/sec. Similar behavior can be found m this situation also The variation of 
has significant influence on zeroth mode and torsional filter mode and its influence on 
exciter mode is small 


Table 6.4: System eigenvalues for different values of a;„ 


(jjn = 2 rad/sec 

Un = 7 rad/sec 

Un = 13 rad/sec 

comments 

-1 38 ± j298 16 
-0 24 ± J202.82 
-0 89 ± jl59 67 
-0 03 ± jl28 93 

-0 08 ± j98 26 

0 01 ± j9 94 
-21.04 ± J20.78 
-1.14 ±jl.61 

-1 38 ± j298 16 
-0 24 ± j202 82 
-0 89 ± jlo9 67 
-0 03 ± jl28 93 
-0 08 ± j98 149 
-2.14 ± j9 18 
-22 32± J20.73 
-2 29 ±j5 25 

-1 38 ± J298.16 
-0 24 ± j202 82 
-0 89 ± jl59 67 
-0 03 ± jl28 93 
-0 08 ± j98 36 
-3 55 ±jl0 38 
-23 29±j20 46 

0 12 ±j7 18 

Mode 5 

Mode 4 

Mode 3 

Mode 2 

Mode 1 

Mode 0 

Exciter mode 
Torsional filter mode 


6,5.2 System with Digital Control based TCSC 

The PSS gain value 4 is selected for the study For this value of gam mode 1 and mode 
2 were unstable without torsional filter The study is earned out foUowmg the procedure 
mentioned in the previous case The results are given in Tables 6 6 and 6.7. From the 
results it can be observed that in this case also the torsional filter is preventing the PSS to 


no 




Table 6.5. System eigenvalues for different values of C 


C = 0 1 

r=05 

C = 08 

comments 

-1 38 ± j298 16 

-1 38 ± j298 16 

-1 38 ±j298.16 

Mode 5 

-0 24 ± j202 82 

-0 24 ± j202 82 

-0 24 ± j202 82 

Mode 4 

-0 89 ±jl59 67 

-0 89 ±jl59 67 

-0 89 ± jl59 67 

Mode 3 

-0 03 ±jl28 93 

-0 03 ± jl28 93 

-0 03 ± jl28 93 

Mode 2 

-0 08 ± j98 26 

-0 08 ± j98 149 

-0 08 ± j98 36 

Mode 1 

- 1 51 ±jl0 94 

-2 14 ± j9 18 

0 055 ± j9 08 

Mode 0 

-22 84 ± j20 78 

-22 32±j20 73 

-22 29±j20 86 

Exciter mode 

0 57 ±j4 1 

-2 29 d:j5 25 

-2 59 ±j5 2 

Torsional filter mode 


interact with the torsional modes The analysis also considered the effect of variation of 
filter parameters From the results it can be noticed that the \’anation of filter paxameters 
are affecting the zeroth mode and torsional filter mode and have minimum influence on 
the exciter mode 


Table 6.6: System eigenvalues for different values of 


= 1 rad/ sec 

u;„ = 5 rad/sec 

ujji = 10 rad/sec j comments 

-85 65± j612 37 
-80 60±jl45 31 
-1 38 ± j298 15 
-0 247 ± j203 12 
-0 88 ±jl59 44 
-0 027 ± J128.87 
-0 07 ± j98 06 
0.23 ± j8 93 
-23.72 ± jl9 78 
-0 38 ± jO.92 

-85 65± j612 37 
-80 60± jl45 31 
-1 38 ± j298 15 
-0 247 ± j203 12 
-0 88 ± jl59 44 
-0 027 ± jl28 87 
-0 07 ± j98 06 
-0 00002 ± j8 955 
-23 78± jl9 80 
-1 69 ± j4 65 

-85 64± j612 37 
-80 60± jl45 31 
-1 38 ± j298 15 
-0 247± j203 12 
-0 88 ± jl59 44 
-0 027 ± jl28 87 
-0 07 ± j98 06 
-3 66 ± jlO 23 
-23 96±jl9 87 
011±j8131 

Supersyn Net Mode 
Subsyn Net Mode 
Mode 5 

Mode 4 

Mode 3 

Mode 2 

Mode 1 

Mode 0 

Exciter mode 
Torsional filter mode 


6.6 Conclusions 


In this chapter, the linearized formulation of IEEE first benchmark model with TCSC, I VR 
and PSS for the study of control mteractions is presented. Both constant angle control 
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Table 6.7. System eigenvalues for different values of C 


C = oi 

CD 

II 

C = 09 

comments 

-85 65± j6l2 37 
-80 60± )145 31 
-1 38 ± j298 15 
-0 247 ± j203 12 
-0 88 ± J159 44 
-0 027 ± jl28 87 
-0 07 ± j98 06 
-1 061 ± jlO 45 
-23 88 ±jl9 86 

0 55 ± )7 56 

-85 65± j612 37 
-80 60±jl45 31 
-1 38 ± j298 15 
-0 247 ± j203 12 
-0 88 ± jl59 44 
-0 027 ± jl28 87 
-0 07 ± j98 06 
-3 1 ± j9 26 
-23 91± jl9 85 
-0 08 ± j8 14 

-85 65± j612 37 
-80 60±jl45 31 
-1 38 ± j298 15 
-0 247± j203 12 
-0 88 ± jl59 44 
-0 027 ± jl28 87 
-0 07 ± j98 06 

0 031 ± j8 63 
-23 98±jl9 79 
-7 77 ± j5 024 

Supersj-n Net Mode 
Subsyn Net Mode 
Mode 5 

Mode 4 

Mode 3 

Mode 2 

Mode 1 

Mode 0 

Exciter mode 
Torsional filter mode 


based TCSC and digital control based TCSC have been considered It can be observed 
from the study, that m the case of fixed series compensation, the AVR is dampmg torsional 
mode 2, while at the same time it is undamping the zeroth mode The mtroduction of 
PSS in this case results in dampmg of the zeroth mode, but it is destabilizmg mode 1 
and mode 2 Higher ^•alues of PSS gam has beneficial effect on zeroth mode dampmg m 
the presence of TCSC but at the same time this results in the destabilization of some of 
the stable torsional modes With the torsional filter all the torsional modes are stable, 
but the variation in filtei parameters affect significantly the zeroth mode and the torsional 
filter mode The filter parameters can be selected such that the maximum dampmg can 
be achieved for the zeroth mode without adversely affecting the system stabihty 



Chapter 7 

CONCLUSIONS 


The objective of the work reported m the thesis is to investigate the occurance of SSR 
phenomenon in a senes compensated power system and to examine the effectiveness of the 
controllable devices in mitigating SSR namely (1) NGH dampmg scheme (2) STATCOM 
(3) TCSC 

7.1 Occurance of SSR 

The occurance of SSR has been studied m case of lEEE-FBM. For this purpose develop- 
ment of system model with fixed senes compensated line is presented m Chapter 2 The 
system model developed in linear state space domain and the modular approach adopted 
in developing the system model permits the representation of a subsystem to any desired 
degree of detail 

The eigenvalue analysis results show that the maclune network interaction in the pres- 
ence of fixed series compensation destabilizes the torsional modes This fact is validated 
and the occurance of instability due to SSR is evident m the time domam results obtained 

7.2 Mitigation of SSR 

The effect of NGH dampmg scheme in mitigating SSR due to fixed senes compensation has 
been examined m Chapter 2 It is observed that although the growth of transient torque 
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oscillations due to SSR are arrested with the NGH damping scheme, but the reduction in 
peak torques is only moderate 

Use of STATCOM has been considered in pow'er system to provide dynamic reactive 
support Coupled with its basic function the effectiveness of STATCOM m mitigatmg 
SSR has been im estigated m Chapter 3 For this purpose a supplementar}- control based on 
reactive current has been considered It is seen that the voltage control alone is inadequate 
in damping the torsional modes This situation significantly improves when the control is 
affected using reactive current along with the voltage control Detailed simulation studies 
also validate this observation and demonstrate the STATCOXI capabilitv with reactive 
current control effectively mitigating SSR while simultaneously providing dynamic reactive 
power support It has also been observed that the STATCOM is more effective in dampmg 
torsional oscillations compared to XGH damping scheme 

Another controllable device which is gaining prominence for enhancement of power 
transfer capability on a transmission line is TCSC While performing this task the TCSC 
IS ('xp('tt('<l to a\oid occurance of SSR The study presented in Chapter 4 demonstrate this 
aspect by considering a constant angle control of TCSC which is based on the philosophy 
of maintaining the voltage drop across the compensated line constant The extent to 
which the power transfer over the line increases would depend upon the level of fixed senes 
compensation considered while the vernier control of TCSC would help m decoupling the 
network modes from the mechanical modes 'While the control design is based on eigenvalue 
analysis, its performance in the event of large disturbances has been quite satisfactory 
For the control of TCSC, from the view point of SSR mitigation, a new digital control 
strategy based on state feedback involving, line current and capacitor voltage has been 
proposed in Chapter 5 To investigate the SSR mitigation capability of TCSC, the model 
of TCSC along with the proposed controller is developed in discrete time linearized domam. 
The results of the analysis show that the digital control scheme makes TCSC extremely 
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effective in damping torsional modes It has also been seen that the TCSC contnbutes 
very significanth to the damping of network modes 

The in\estigations on the effectneness of the various controllable devices mentioned 
above did not consider the presence of excitation system The design of a controllable 
device (STATCOM or TCSC) should generally consider the presence of any other control 
present in the sjstem primarily to ensure that the two controls do not have conflicting 
contributions to the system performance 

The interaction between the generator excitation system (AVR and PSS) and TCSC 
has been studied and it is observed that while PSS has favorable impact on the zeroth 
mode (power swing mode), and it has a detrimental effect on the torsional modes This 
impact of PSS conflicts with that of TCSC This situation can overcome with the use of 
torsional filter along with the PSS has shown in Chapter 6 The torsional filter makes the 
PSS insensitive to torsional modes however, the parameters of torsional filter need to be 
carefully chosen to ensure overall system stability 

7.3 Suggestions for future work 

1 The effectn eness of TCSC or STATCOM in mitigating SSR need to be examined 
in multi machine environment This investigation also result in alternate control 
strategies 

2 Application of the generator rotor speed synthesized from locally measurable quan- 
tities, as supplementary signal m mitigating SSR can be examined 

3 The frequency scan method provides information regarding possible problems with 
torsional interaction and transient torques The study can be made with the fre- 
quency scanning technique to further examine the effec iveness TCSC m mitigating 
the SSR 
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Appendix A 

Generator System Matrices 


In this appendix generator system matrices are given 


A.l Synchronous machine model matrices 


The Synchronous machine model matrices are given below 


[Asi] 


Ws 

0 


1 

T' 

90 





-id. 




d J 


Bjii 


^B^gq —OJB'^q 

0 0 

0 0 


[B,i 3] 


Bju] = 


0 
0 
0 

L i 


—ubcosS 

—cjBsmS 


0 

0 


UBStnd 

—ubcosS 

0 

0 
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For the six mass mechanical system non-zero entries of the matrix and are 

given below 

Am IS a 12 X 12 matrix 


.4r,»(l,l) 
-4„^(1,2) 
-4.„(1,7) 
-4m(2,l) 
^4^(2, 2) 
A.. (2, 3) 
-4m (2, 7) 
-4m (2, 8) 

-4m (3, 2) 
.4m (3 3) 
-4m (3, 4) 
-4m(3,8): 
-4m (3, 9): 
■4.m(4, 3) : 
-4m (4, 4) : 
-4m (4,5) = 
Am(4,5): 
■4m (4, 10) 
4m(5,4) : 
4m(5,5) = 
4m (5, 6) = 


{D.+D. 






2Hc 

1 

2H, 

2Ha 


{Dq^Dqe+DgLPB) 


2Hg 

— - ^9LPB 

— 2Hg 

L 

2Hg 

2Hg 

= ^9lPB 
2Hcpb 

= (^IPB'^DqLPB + ^LPBLPa) 

2Hipb 

— pblpa) 

^ 2Hipb 

_ I 

2Hipb 

— I 


2Hipb 

{Dlpbipa) 

2Hipa 

{Dl PA^DlPALPB+I^tPAIp) 

2HtPA 

Olpmp) 

2Hlpa 

I 


2Hipa 

- 1 
2Hipa 

Dlpaip) 

2HiP 


(Djp’j-DLPAIPi’DiPHP) 

2HiP 


Diphp) 

2Hip 
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10 ) = - 


1 

2H[p 


^4^(5, 11) 


2Hjp 


^4r„(6, 5) 


2Hhp 


4^(6, 6) = -^-^^i^=£^Hp±£EPl 


-4m(6, 11) — — 2 Hpp 


- 4 ^( 7 , 1 ) = -Kg, 

Ami 7 r 2 ) = Kg, 

-4m(8 2) = —hgipQ 

-4t7i(8, 3) = hgiPB 
-4„i(9,3) = —Klpblp i 
.4m (9, 4) = Kipblpa 
.4m(10, 4) = -Kip^jp 
-4m(10,5) = Kipj^ip 
-4m(ll,5) = —Kipbp 
Ara{lh%)=-KiPHP 
.4m(12, 2)=Ub 

[Bm] IS a 12.Y1 matrix with non-zero entry 

^m(2 1) = ~ 2 W^ 

[Cm] IS a 2.Y12 matrix with non-zero entry 

Cm(l,12) = l 

Cm (2, 2) = 1 


The matrix [Cgi] is defined as 


[Cgi] = [Ccii] AXsi 4- [C^Glm] AXm 
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wheie 


o 

o 

09 

II 

cos5 

sinS 

sin5 

cosd 

Sind 

cosd 

ctls 


1 

. ^'d 


^0 



For a SIX mass s}steni 

[C(71)7i] is a 2 jY 12 matrix with non-zero entries 
C7Glm(l, 12) = iQo 
C'Glm(2, 12) = -iDo 


The s\nchronous machine output matrices are given below 


[Cg] = 


Cm 

CgiAg 


and [ Dg ] 


0 

CgiBg 


In Chapter 5 the following synchronous machine output equation is used to get the overall 
system model 


Aye = [Cg] Axg 


where Ay a = [Aio Aiq] and [Cg] = [Cgi] 


The interface matrix [F] for the lEEE-FBM is given by 


[F] = 


Ri Xi Ti 0 1 0 

-AT 0 Li 0 1 

1 0 0 0 0 0 

0 1 0 0 0 0 
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Appendix B 

IEEE First bench mark model data 


The data for the Synchronous machine is adapted from IEEE first benchmark model [1] 
The data is given on 892 4 MVA base 
(i) Generator data 
/ = 60 Hz, Ra = 0 

Xrf = 1 79 per unit , = 0 169 per unit, = 4 3 sec x, = 1 71 per umt. = 0 228 per 

unit, = 0 85 sec. 


Table B.l; Mechamcal system data 


Mass 

Shaft 

H (second) 

K (per unit) 

HP 

HP-IP 

0 092897 

7277 

IP 

IP-LPA 

0 155589 

13168 

LPA 

LPA-LPB 

0 85867 

19618 

LPB 

LPB-GEN 

0 884215 

26713 

GEN 

GEN-EXC 

0 868495 

1064 

EXC 


0 0342165 



(ii) Transmission line data 

Rl = 0 02 per unit , Xl = 0 5 per unit 
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(ill) Transformer reactance 

Xr = 02 per unit 

(ill) Operating point data 

Vg = l per unit P = 0 9 per unit , at 0 9 p flag 
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Appendix C 


Synchronous machine initial 
condition calculation 


The vectoi diagram for the generator is shown in Fig C 1 The generator terminal voltage 
I'ZS and current r,Ze‘ is known m phtaor form from the initial conditions of the network 
variables and power flow through the transmission lines The rotor angle 6 and the initial 
values of £?', E'^ and -tpg are calculated as follows 


EqZS = VgE9 -i- IgEd' {JR^ -J- jxq) 
Vgg = VgCOS{S “ 

Vgd = -VgSm{6 - 9) 
iq = IgCos{5 — 9 ') 

Id = —IgSin{5 — 9 ') 

Eq — Vgq Xd}d d" Ra^g 
Ed — Vgd d” Ra^d 

= x’dld d- Eg 
'tpg = Xqlq - E'd, 
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where and are the dg components of the 
corresponding current quantities 


generator ’toltage and and 



q axis 


Q 


axis 


Figure C 1: Synchronous machme phasor diagram 



Appendix D 

STATCOM Controller and 
Interconnection Matrices 


Thi!. appendLX deals wUh the STATCOM centroUer matrices interconnection matrices 
initial condition calculation and STATCOM data 


D.l STATCOM Controller Matrices 


The matrices related to STATCOM controller are given below 


where 


[Asc] = 


0 0 0 ' 
-if/ 0 0 0 

0 0 -ft 0 

-K^K2 K2 -K2 0 , 


[Bsci] = 


_K^ KoV,o 
V,T,n 

0 0 

ik_ 

VsTm ViTm 

0 0 


[Bsc 2 ] = 


iil 2 _ 

V,T„ 


Ejisi 

Tm 


Tm 

0 


_ Ej 2<^ ■ 

VsTm Tm, 

0 

Si. 

Tm 

0 


^sQ ‘iRYsP 

V, Vi 
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and 


C 3 


= !f£ _ ’’'R^sq 




K2 


[^sc] = [~K,K, K, -K, 1 ] 
[Cci] = f 1 0 0 0 0 0 0 ^ 

[o 1 0 0 0 0 0 

[^sa] = [l^Vouf all columns of As, 2 "^ 


row all columns of A5] 


[C 5 ] = 


C51 

C52 


[Dsi] IS null matruc of size (2x2) 

[0^2] = [l“ r<TO all columns of Bj, 2”' row all columns of Bs] 


[D 5 | = 


Dsi 

D52 


D.2 Details of interconnection matrices 


The [F] matrix describes the interconnection between inputs and outputs of vanous sub- 


systems, is given by 


m= 


Til Fi2 
T21 F22 
T31 F32 


Ti 3 

T23 

F33 


The elements of the [F] matrnx are mentioned below 


[Fii] = 


[F 12] = 


Ri Xl Li 0 

—Xl Rl 0 Lz 

Ri Xl Li 0 
-Xi Rl 0 Ll 


[F13] = 


1 0 
0 1 


[F21] 


Rl Xl Ll 0 
—Xl Rl 0 Ll 
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[^ 22 ] 


Rl Xl Li 0 
—Xl Rl 0 Li 


[F23] 


1 0 ' 
0 1 


[F31] 


10 0 0 
0 10 0 


[F32] = 


-10 0 0 
0-100 


[F33] 


0 0 
0 0 


D.3 STATCOM Initial condition calculation 


The procedure for caluculating the imtial conditions is mentioned below From load flow 
studies, Pj, V", and Qs at STATCOM bus are obtained 
Bs which is the susceptance of the STATCOM is given by 



The current h through the STATCOM is given by 


^5 — jRs'^S 


The DC \oltage Vdc is given by 


Vs + 

Vdc = Z 


and the angle /? is given by 


^ = Xvs- X {vs - ^s{R + jxi)) 
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Appendix E 

TCSC Controller and 
Interconnection Matrices 


This appendix deals with the network matrices, demotion of contmuous time TCSC model 
given in Chapter 5 and interconnection matnces 


E.l MATRICES FOR TCSC CONTROL 


The network matrices are given below 


B 


N 




0 

LJg 

0 

0 


-UJQ 0 

0 0 

0 0 

0 (jJb 


ojb Xcf 
0 

XtcSC 

0 


0 

ub Xcf 
0 


0 

0 

— u^B 
0 


~'^TCSCD -“-l UilU n, u u I 

wb Xtcsc -■'^tcsco Xtcsc A: 0 0 J 

dXrnsce ff(cr) 


0 0 0 

0 0 0 

^TcscD '^B Xfosc ^00 

» f.ln fl 0 


where 


k 


7.2 

A/O 


cosa , „ , s . 

{Xc + ki) + 


IT 


TTXo 


cosa , ^ f k2<J 

kostna tan —— 

2 V 2 


cos 


(I) 


sec? (^) 


XTCSC,„io)=X,- ^'^'^^^ lXc + hl + ^cos^{^) [tan(l) 


k2a 
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D 


V 


k, = 





ko — 


ks — Xxcsceff{(y) 

1 0 1 o' 
0 10 1 
0 0 10 
0 0 0 1 


'N 


Rl ws Ll 0 Li 0 
-wb Li Ri 0 0 Ll 

0 0 0 0 0 

0 0 0 0 0 


E.1.1 Derivation of TCSC Voltage 


The TCSC model is developed based on dynamic phasor approach In this approach a 
waveform t(.) can be represented on the interval {t — T,t) usmg a Fourier series of the 
form 

x { t ) = Re [EXfc(t)e^-'»^] 

wheie k > 0, Re denotes the real part, ws = ^ Xk are the Fourier 

coefficients, which are referred as phasors These Fourier coefficients are functions of time 
since the iiitorval under consideration slides as a function of time The k-th. coefficient (or 
k — phasor) at time t is determined by the following averaging operation 

X,(t) = I r x(T)e-”^’'dT =< X >, (t) 

Jt-T 

where c = 1 for A: = 0 and c = 2 for A: > 0 The den^•atlve of the A;-th Fourier coefficient 
is given by the following expression 

at at 
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The TCSC capacitor voltage and parallel inductor current are denoted by vjcsc and Zp 

respectively The parallel inductor is denoted by Lp The line current is denoted by zg 
T he state space model is given by 


C 


dip 


where 9 is a switching function that denotes the thynsior status ? = 1 when one thsTistor 
IS on, and 5' = 0 when both thyristors are off The Ime current is assumed sinus-jidad over 
any interval of length T = ^ The state space model is obtained by applymg averaging 
operation to the above equation resulting in 


.dv. 


CTCSC^^ = ^o- zpx - jujbCtcscv. 


dt 
^dt 
dt 


TCSCl 


~ir =< <l'^Tcsc >1 -y^^BLpipi 


The quasi steady state model of TCSC is obtained with the following assumption 

'^tCSCI 


tpi 


where 


J^BLeffia-} 

1 


UJb^{CtcSC - Crcsceffia)) 

The quasi steady state model of the TCSC is obtained using the above approximation as 


C'TCSC — = iQ — (jCBbCtCSC d r )‘^TCSCi ~ J^BCpcSCeffia ^tcsci 

In the above equation the TCSC voltage and line current are fundamental phasor quantities 
1 e , Ifc = 1 The above equation is derived in synchronously rotating reference frame and 
can be written as 

where, a = cr^ef + 2 arg {-Pg'^tcsc) 
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E.2 Details of Interconnection matrices 

The [F] matxix which describes the interconnection between inputs and outputs of the 
various subsystems is given by 


where, for the study system F 


11 



Fii 

Fi2 

Fi3 

[Fl = 

F21 

F22 

F23 


. F31 

F32 

F33 

= Fi3 

= 

f'22 = 

= F33 = 0 



■ 1 

0 0 

0 ■ 



0 

1 0 

0 

[F2,I = 

= 

0 

0 0 

0 



0 

0 1 

0 



_ 0 

0 0 

1 _ 

[Pl2l ^ 


■ 1 

0 0 

0 ' 


_ 0 

1 0 

0 




■ 0 ■ 





0 


[F23 

] = 

1 





0 





_ 0 _ 




■ 0 

0 0 

0 ' 

[Fail = 


0 

0 0 

0 

- 

1 

0 0 

0 



_ 0 

1 0 

0 _ 



■ 0 

0 1 

0 ' 

[F32] = 


0 

0 0 

1 


0 

0 0 

0 



0 

0 0 

0 
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Appendix F 

Excitation System and PSS Matrices 


In this appendix the matrices associated with chapter 6 and data of exciter and PSS are 
mentioned 


F.l Exciter and PSS Matrices 


The matrices related to excitation system smd PSS are given below 


Agi — 


i_ 

0 

0 

0 


Bell — 


Tr Tr ^2 ■ 

Kg^Cl 1 

-Cl C4 

K.C, 

0 -Cl 

KsC^ 

0 0 

1 

Tu, J 


■ ifcqKg - 



KsCs 



KgC^ 



1 

L J 



Bel2 = 


KRVgP 

TRVg 

0 

0 

0 


KrV(!Q 

TRVg 

0 

0 

0 


Eel = 


M 1 

Tr 

0 

0 

0 


Cel = [ 1 0 0 0 ] 
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where C. = X, C. = c, = 1 - 1, c, = c.Ca, C. = CrC.q 
time constant 


3 and Tw is the washout 


F.2 Exciter PSS and Torsional filter matrices 

The matrices related to excitation system, PSS and torsional filter 


Agi — 


__i_ 

Tr 

0 

0 

0 

0 

0 






are given below 


0 

0 

0 


0 

1 

0 

0 

0 


0 

0 

-Cl 

0 

0 


0 

0 

4C2 

C4 

-Cl 

0 


0 

0 

^IClKs 

CiKs 


Bell = 


0 

0 

K.C, 

KsC 4 

1 


B 


el2 


^RVgD KrV^q 1 


TRVg 

0 

0 

0 

0 

0 


TRVg 

0 

0 

0 

0 

0 


Eel = 


r Aa 

Tr 
0 
0 
0 
0 
0 


Cel = [ 1 0 0 0 0 0 ] 
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where C. = X, Ce = c, 
the wahout time constant 


1 - H 
r2> 


<^4 = C1C3, C5 = Cl C, C3 and Ti 


ly IS 


F .3 Exciter and PSS Data 


The excitation system and PSS data is given below 
Exciter Data 

Kr = 200, = 0 025s, Tw = 10 

PSS Data 

Ti = 0 048s, T 2 = 0 032 
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